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ABSTRACT 

This report summarizes research done during the past twelve months 
under NASA Grant NSG-1116, entitled A Theoretical Study of Hexero junction 
and Graded Band Gap Type Solar Cells. The purpose of this work is to 
design a computer program for the analysis of variable composition solar 
cells and apply that program to several proposed solar cell structures 
using appropriate semiconductor materials. As -discussed in this report, 
the program has been completed and has been used to study devices made 
of A1 Ga, As, Ga In- As and GaAs n P with the intention of determining 
the material compositions, and device dimensions that produce high 
conversion efficiency. 

The computer program has been designed to simulate solar cells made 
of a, ternary alloy (such as A1 Ga. As) of two binary semiconductors 
with an arbitrary composition profile, and an abrupt or Gaussian doping 
profile of polarity n-on-p or p-on-n with arbitrary doping levels. An 
antireflection layer of SiO of any reasonable thickness can be 'Specified 
and either AMO or AM2 solar irradiance conditions can be selected. Once 
the device structure is specified, the program numerically solves a com- 
plete set of differential equations 'and calculates electrostatic 
potential, quasi-Fermi levels, carrier concentrations and current 
densities , total current density and efficiency as functions of terminal 
voltage and position within the cell. These- results are then recorded 
by computer in tabulated or plotted form for interpretation by the user. 

The report includes background discussions concerning variable com- 
position solar cell phenomena and the theoretical techniques used to 


' Preceding page blank 



model device behavior. The results of the computer analysis of various 
proposed cell structures are then presented. It was found that inter- 
face recombination due to lattice mismatch at abrupt AlAs-GaAs hetero- 
junctions produced an insignificanr reduction in peak efficiency while 
recombination at GaP-GaAs heterojunctions produced far more severe 
losses. However, it was subsequently found that these losses could be 
largely eliminated by placing rhe p-n junction slightly below the GaP- 
GaAs hetero junction (by about 0.04 yM) and by grading the composition 
change over a small distance (also about 0.04 yM). 

According to computer simulation, even a very thin window layer 
(0.1 yM) of AlAs or GaP should virtually eliminate the surface recombina- 
tion losses known to limit GaAs homojunction solar cells. In fact, to 
minimize generation-recombination loss in the window and maximize 
overall efficiency, the wide bandgap layer should be as thin as tech- 
nologically feasible (at least down to 0.1 yM) . Finally, the problem 
of optimizing the composition of the substrate in AlAs-on-Al Ga. As, 

5C -LX 

AlAs-on-In Ga As and GaP-on-GaAs P structures is explored by 

X X — X X” X X 

simulating cells with various substrate compositions. 

These computer studies are intended to provide a more accurate 
estimate of the potential performance of variable composition solar cells 
than is possible with simpler, non-computer models. The results of these 
simulations should serve as a guide for the design of high efficiency 
experimental solar cells. 
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1. INTRODUCTION 


1.1 Incentive for Computer* Simulation 

Recent developments in the world's energy market have underscored 
the need for development of alternative energy sources to reduce our 
national dependence on fossil fuels and provide for the Inevitable 
depletion of such limited resources. Although solar energy cannot be 
expected to satisfy a significant fraction of the world's energy demands 
for several decades, its unlimited nature will make it progressively 
more competitive as scarcity forces the cost of existing sources higher. 

Given the abundant supply of solar energy and the growing demand 
for electric power, the ultimate answer to the energy problem would 
appear to be the solar cell, which directly converts sunlight to elec- 
tricity. However , the high cost per watt (currently about , $15 . 00 ) 
associated with solar cells has limited their use to a few specialized 
applications Cl, 60 j.-> -Reduction of ‘this price depends -on the 1 reduction of 
cell fabrication costs, and an increase in conversion efficiency, which 
typically ranges from 12% to 16% depending on cell structure and 
illumination conditions . 

While most solar cells are presently made of a single semiconduc- 
tor (usually silicon), some of the most efficient cells have been built 
by allowing the material composition to vary with depth through the 
device so that various layers can have independently controlled optical 
and electrical characteristics. Such a cell is diagramed in Figure 1.1, 
where the cell substrate consists of GaAs, and the surface layer is a 
linearly graded alloy of A1 Ga. As. Preliminary calculations predict 
a conversion efficiency ’of approximately 17% under air mass zero (AMO) 
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Figure 1.1 General composition profile and structure for a 

variable composition A1 Ga As ^olar cell. 
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conditions for* such a device compared with 9% for a good homojunction 
GaAs cell [49]. Although it is to be expected that variable composition 
cells will carry higher fabrication costs per unit area, the improvement 
in efficiency could make them attractive in applications where other 
costs dominate the complete system price, such as those having concen- 
trator lenses and solar tracking mechanisms. Variable composition 
A1 Ga As cells potentially have a high conversion efficiency because 
they can prevent the severe surface recombination losses associated 
with GaAs homojunction cells while retaining the superior match to the 
solar spectrum which constitutes the main advantage of GaAs cells over 
Si cells. 

Since the design of variable composition solar cells permits xhe 
use of a wide variety of semiconductor alloys and composition profiles , 
it would be useful to have a computer model that could accurately pre- 
dict the performance of proposed cell structures before they are 
actually built. For example, it is not obvious whether maxerial 
composition should be varied gradually wirh depth to create a graded 
layer (as in Figure 1,1) or abruptly to create a hetero junction. Like- 
wise, the determination of the optimum placement of the p-n junction in 
relation to the other layers is not an easy problem. In an effort to 
answer these, and other questions, an existing homojunction solar cell 
computer analysis program has been modified to permit the modeling of 
variable composition solar cell structures. The computer model serves 

i / 

as an intermediate step between first order approximations of device 
behavior and the physical construction of a working model. The computer 
program can include the many detailed second-order phenomena that must 



be excluded from first order calculations,. In addition, the structure 
and specifications of a computer model can be manipulated with a flexi- 
bility that is impractical during the testing of a physical. model. Thus, 
computer analysis can point out the critical design factors and help 
determine the most efficient and least expensive structure before it is 
actually built. 

The present work describes the ideas involved in a computer program 
development and the performance it predicts for various proposed cell 
structures. Also included is a discussion of the relative importance of 
various design factors in obtaining, maximum efficiency based on computer 
predictions and the limitations of solid .state technology. 

1.2 Literature Review 

\ 

The development of analysis' techniques for variable , composition 
solar cells can be traced back to attempts by several workers to explain 
early -silicon homo junction solar cell .behavior . In 1955', Prince [20] 
approximated the J-V characteristic of such cells by shifting the well 
known Shockley diode curve by the amount of short circuit 'current that 
could be collected from optically generated carriers. An estimate of 
internal resistance losses was also included in the model. In 1956, 
work by Loferski [21] gave the first indication that silicon. did not 
provide the optimum optical 'match to the .solar spectrum under AMO conr 
ditions. In fact, in the same year, Jenny et al. [57] reported an 
experimental GaAs solar cell efficiency of 6.5% which was comparable 
to the values obtained from Si and CdS solar cells at that time. 

However, due to the dominance of silicon technology in fabricating other 
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solid state devices, most 'solar cell experimentation continued to con- 

i 

sisx of variations on basic silicon structures. 

The first theoretical examinations of graded bandgap materials in 
solar cells came in 1957 when , Tauc [3] estimated the photovoltage due 
to bandgap variation. In 1960, Wolf [4] analyzed the behavior of a 
graded gap region in combination with a p-n junction, and later, in 
1962, Em t age [53 predicted that such a configuration could produce 
conversion efficiencies -as high as 35%. At about the same time, 

Anderson [55] developed a basic model for heterojunctions which 
considered current flow by the -injection of carriers over band edge 
barriers j and later, in *1964, the abrupt heterojunction solar cell was 
theoretically analyzed by Perlman. and Feucht [56]. 

The suggestion by. Kroeraer in -1963 [9] that the use of heterojunc- 
tions (interfaces between dissimilar materials) could improve the 
performance of injection lasers was followed by escalated laser research 
using position-dependent, ternary III-V .alloys 'for light and carrier 
confinement. The resulting need for experimental data concerning the 
material properties and band structures of III-V alloys was satisfied 
by many workers, notably Cohen and Bergstresser [7] in 1966 and by 
Thompson and Woolley [8] in 1967, Such data would prove vital to later 
research on variable composition solar cells. 

By the late 1960’s , various band diagram models for semiconductor 
heterojunctions , had been proposed and in 1971 Marfaing and Chevallier 
[10] published a theoretical analysis of photovoltaic effects in 
variable composition devices that allowed -for position dependent band- 
gap, carrier mobilities, effective masses, lifetimes and doping. 
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Their theory was accompanied by their experimental studies -of Cd Hg Te 

X -L ““X 

graded composition crystals. In 1973, Wemac and Rediker [11] explained 
hetero junction laser behavior by ’Using a model 'that spread the composi- 
tion variations over- a few hundred angstroms, resulting in improved 
understanding of carrier behavior in all' variable -composition devices. 

At about the same time, 'Hovel and 'Woodall [12] achieved a solar cell 
conversion efficiency of 14.7% tinder _ AMO /conditions by placing an 

A1 Ga- As window layer over a GaAs substrate. This represented a sub- 
X X*“X 

stantial improvement over the 9.3% efficiency reported- in 1972 by 
Tsaur et al. [22] for GaAs cells without 'a wide gap window layer and. 
operating under similar conditions. 

The computer techniques used in the present research stem from 
work by several predecessors. In 1970, Graham and Hauser [13] 
developed a computer analysis program for p-n junction diodes by 
building on techniques Used by earlier workers, notably Bellman and 
Kalalba [14,15,-16], DeMari [17] and Lee [18]. Four years later, Dunbar 
and Hauser [19] modified the program so that 'it could simulate silicon 
homo junction solar cells, The essential purpose of these computer 
studies was the analysis of solar cell behavior without resorting to 
first order approximations such as those necessary for the solution 
of nonvcomputer models. The present -Work represents -an extension of 
this idea into the more. complex field of variable composition solar 
cell analysis. 
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1.3 Computer Analysis Techniques 

The compuxer program used to model solar cell behavior is based on 
the familiar set of nonlinear differential equations normally applied xo 
p-n junction diodes. These consist of the transport and continuity 
equations for electrons and holes, along with Poisson's equation. The 
influence of light. on the diode equations appears as an independently 
calculated optical generation rate, G^, inserted in the continuity 
equations. The material parameters, such as mobility, lifetime and 
bandgap, are provided as functions of composition and/or doping (or 
position if the composition and doping profiles are given). An attempt 
has been made to model all parameters as accurately as possible by a 
combination of theoretical expressions and numerical data. Boundary 
conditions concerning surface recombination velocity and terminal 
voltage are then applied to the complete set of equations. 

In order to solve this system of nonlinear differential equations, 
xhe program uses a technique known as quasi-linearization used 
previously for diode analysis and later used in homojunction solar cell 
simulation [13,19]. Briefly, this technique replaces the non-linear 
equations with linear approximations derived from truncation of the 
expansion of nonlinear terms. The linear equations are then solved 
iteratively to converge on the solution to the original nonlinear 
equations. The results consist of -the printout and plotting of electro- 
static potential, and current and carrier densities for elecxrons and 
holes as functions of position within the solar cell. The most signi- 
ficant results of these analyses are discussed in later chapters of 
this report. 
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1.4 Organization of the Report 
The following chapters deal with general background discussion, 
specific computer techniques, and significant analysis results obxained 
from computer simulations. 'Chapter 2 begins with an explanation of 
homojuncxion solar cell operation and leads into ideas which explain 
variable composition solar\ cell behavior . Chapter 3 presents the 
basic homojunction device equations and develops the modified set of 
equations used to study variable composition cells. The detailed 
methods used to incorporate material parameter data into the program 
are illustrated in Chapter 4 because the assumptions inherent in 
those methods can be crucial to proper interpretation of the analysis 
results. Chapter 5 describes the semiconductor materials that were 
modeled along with the reasons for their selection. In Chapter 6, the 
most significant of the results obtained from the computer simulation 
of these structures are presented and interpreted in light of the 
analysis techniques and the limitations of semiconductor technology. 
Finally, Chapter 7 summarizes the important ideas associated with this 
work and recommends areas of future research. 
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2. PRINCIPLES OF SOLAR CELL OPERATION 

2 . 1 Homojunction Solar Cell Behavior- 

Figure 2,1 shows a typical n-on-p silicon solar cell having an SiO 
antireflectiop layer at the front. Solar photons pass into the n 
layer with a transmission efficiency enhanced by the SiO layer. As 
they pass through the silicon crystal, those photons having energy 
greater than the bandgap of silicon excite electrons across the gap 
generating excess electron-hole pairs . The built-in field set up by the 
doping gradient around the p-n junction separates these excess holes and 
electrons creating a potential across the device terminals that can be 
used to force a current through an external load. However, since none 
of these processes is 100% efficient, most of the solar energy striking 
the cell surface (typically more than 80%) fails to be converted to 
electrical energy delivered to the external load [19] . 

First of all, since the solar spectrum contains a range of photon 
energies (see Figure 2.2), no antireflection layer can insure trans- 
mission of all photons into the cell. For the same reason, some photon 
energies will be smaller than -the bandgap, and pass through the cell 
without generating a significant number of electron-hole pairs. Other 
photons, having energies greater than the bandgap, will generate 
carriers but the energy in excess of that necessary to bridge the gap 
will eventually serve to heat up the crystal. Since these optical 
losses depend mainly on the semiconductor bandgap, their manipulation 
is severely limited If the cell consists entirely of a single semi- 
conductor and therefore a constant bandgap. As will be discussed in 
later sections, the introduction of position-dependent composition 



SiO 


Photons '/•: 





Figure 2.1 Typical n-on-p silicon solar cell. 
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provides an additional degree of freedom, in minimizing these optical, 
losses which often .amount to over 50% of the .incident solar endfgy. 

The second set of loss mechanisms is due to a failure. to collect 
all the excess carriers that are generated within 'the device * The , 
energy captured by generated carriers -can be lost as electrons and 
holes recombine in the junction depletion regio.n 9 in the bulk,- or at 
imperfections at the front surface* In a homo junction cell, the 
collection of excess minority carriers, is due to the built-in electric 
f ield , created by the ;doping gradient around the p~n junction. The 
efficiency with -which the field sweeps these .carriers .across the. 
junction depends on the distance between the junction and their point 
of generation and upon their diffusion length, which is in general a 
function Of position. 

From the well known transport equations for carrier behavior in 

L \ 

diodes, the built-in field in the quasi-neutral -regions (away from the 
p-n junction depletion region) at equilibrium is approximately: 


. w d <* n y 

q dx 


(n material, R 


D 


n. ) 


E = 


^ d(£n 
q dX 



material, 


» n. ) 

i 


Thus, doping gradients -in- the bulk can be .used to help Sweep minority 
carriers to the p-n junction depletion region where the built-in field 
is even greater. Although it would be desirable to have the region of 
ffiaxiffitim optical generation .rate in "or near the region Of maximum field 
(around the. p-n junction), the achievement ,of SUch a situation in a 
honJojuncnipn cell requires a very thin surface layer. But, a thin 
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surface layer not only presents a high series sheet resistance, but also 
places the front surface (with its short carrier lifetime) near the 
junction region as well. In fact, as -long as a single material is used 
throughout the cell, it is impossible to separate the region of maximum 
generation rate from the region of maximum recombination simply because 
the generation rate of a homo junction cell (see Figure 2.3) decreases 
monotonically with depth beneath the front surface. As will be shown 
later, the introduction of position-dependent composition allows -the 
placement of peak generation well below the front surface. 

The final group of loss mechanisms results from the electrical 
behavior of carriers after they are collected. First of all, as 
minority carriers cross the p-n junction they become majority carriers 
and unless short circuit condixions exist, they reduce the built-in 
junction potential and forward bias the diode. The resulting forward 
current constitutes a loss that is determined by the dark J-V charac- 
teristics of the cell. Such an effect can be modeled to first- order by 
the circuit shown in Figure 2.4a where the diode element characterizes 
solar cell behavior in the dark, and the current density source equals 
the short circuit current delivered by the cell under specified 
illumination conditions . This relationship can be written as 

1 t 

J d = J SC “ J S C2,l) 

The open circuit voltage, V QC > occurs where the current generated by 

solar illumination exactly offsets the 'current lost as forward bias 

1 

current in the diode. 


V 

oc 


mkT 
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( 2 . 2 ) 



GENERATION RATE (#/CC/SEC) 





16 


Because of the nonlinear diode characteristics, the maximum power obtain- 
able from any solar cell is determined not only by the short circuit 
current and open circuit voltage, but also by the sharpness of the 
"knee" of the J-V curve as measured by the curve factor: 


C 


f 



(2.3) 


sc oc 

The curve factor is a complex function of the dark characteristics of 
the solar cell diode. 

This simple model leads to an expression for maximum power density 
(excluding series resistance losses) 


P . = (J C_) (— £n-~) 
peak sc f q J g 


(2.4) 


If the reverse saturation current density, J , is due predominantly 

s 

to diffusion in the p-type base (an n -on-p structure will be assumed 
for discussion purposes) then the factor, m, in the diode equation 
approaches 1, and J can be approximated as 


J = 
s 


q/T n.' 
^ n i 

/r N 


(2.5) 


where D is the electron diffusion coefficient, n. is intrinsic carrier 
n l 

density, t is excess carrier lifetime in the base, and is the doping 

density in the base. If J g is controlled mainly by depletion region 

recombination, then>m approaches 2 and J can.be expressed approximately. 

s 

by 


J = 
s 


q W,n. 
d ,i 


x + T 
po no 


where is the width of the depletion region. 


( 2 . 6 ) 
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A particular 1 solar cell 'may -have, characteristics that 'lie somewhere 
between these .two situations but. in any case it is apparent that the . 
task of maximizing peak power output is' full of tradeoffs. For , examples 
from .Equation (2,4), it appears that in order to -maximize peak .power 
density, should be minimized, possibly by minimizing and B^ (if 
Equation (2.5) applies). But, this would .reduce J sc since 'a smaller B^ 
means a lower carrier collection .efficiency, and a smaller would 
generally be achieved. by choosing a material with a .larger bandgap, in 
turn making the crystal Transparent to a larger share of, the solar 
spectrum. Similarly, although Equation (2.4) apparently indicates that 
? peak directl y Proportional to m, -C f decreases, with increasing m, so 
that the overall influence < of m on peak power is -not obvious . from . 
Equation (2.4). In short, the influence of .several design factors. on 
solar cell performance cannot be easily predicted. By computer .simula- 
tion using accurate modeling techniques and reliable data, the dominant 
factors controlling solar cell behavior .can be, clarified . 

A final loss_ factor is associated with the layout' of contacts on 
the front surface. As can be seen from Figure 2.5, their 'design involves 1 
a compromise between minimizing series sheet resistance loss due to 
lateral current flow, and maximizing the surface exposure to .light. 

The computer, program ~ does .not account for these resistive and optical 
losses because irhey can be reasonably estimated for any cell once 'a 
contact geometry ,is -.specified, and because the aim of this -research is 
to determine the relative influence of various .material parameters 'that 
control the 'internal operation of the cell; It should be noted,- however, 
that internal bulk resistance loss is accounted .for inherently in the 



18 


photons 



by typical finger contact grid. 
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Table 2.1 Summary of Solar Sell loss mechanisms. 

Optical Losses 

1. Surface reflection 

a) From metal coverage 

b) From exposed semiconductor 

2. Photon energy in excess of bandgap in excess of 

that necessary to generate electron-hole pairs) 

3. Failure to absorb all photons because t 

a) Photon energy less than the bandgap 

b) Cell too thin to absorb essentially all photons 
of some wavelengths 

Collection Losses 

1. Surface recombination (front and bach) 

2 „ Bulk recombination 
Electrical losses : 

1. Internal diode bias current 

2. Leakage current at imperfect cell edges 

3. Series resistance voltage drops due tos 

a) Bulk resistance 

b) Contact resistance seen by lateral current flow 
to contact fingers 



program by the device equations discussed in Chapter 3. Table 2.1 
summarizes the important loss mechanisms that dominate solar cell 
behavior . 
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2.2 Variable Composition Solar Cell Concepts 
2.2.1 Introduction 

The phenomena that are responsible for variable composition solar 
cell behavior are essentially the same' as N those found in homo junction 
solar cells. The major distinctions between the two stem mainly from 
the wide variety of variable composition structures which are available 
and from the increased complexity of the relationships between energy 
conversion and loss processes that control performance. However, this 
complexity brings with it an additional degree. of freedom in- controlling 
losses and maximizing conversion efficiency. 

For discussion purposes, cells made of AlAs alloyed with GaAs in 
varying proportions will be examined in this section. As shown in 
Figure 2.6 the alloy proportions can vary gradually to form a -graded 
layer, can be changed abruptly ,to form a heterojunction or any number 
of other variations. The composition at the front surface, and in the 
bulk can be arbitrarily specified, and the p~n junction may or may not 
be placed within a graded region. 

The relative merits of" these, and other options will be discussed 
conceptually in the following sections and the computer simulation 
results to be presented in later .chapters provide a basis for evalua- . 
tion of these concepts. 
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2.2.2 Generation Rate 

The most ' outsxanding impact of position dependent composition is on 
the optical generation rate profile. The profiles of a homojunction 
GaAs cell, a linearly graded A1 Ga As cell (see Figure 2.6a), and an 
AlAs-GaAs hetero junction cell (see Figure 2.6b) are shown in Figures 
2.7, 2.8 and 2.9 respectively. . The calculations used to determine 
these curves 'are explained in Chapter, 4. 

The profile of a 1 homo junction cell consists of the sum of , the 
individual generation rates -due to absorption at each wavelength in the 
solar spectrum. Since the absorption 'coefficient for a given wavelength 
is constant if material composition is constant, the generation rate for 
each wavelength decays exponentially, and the ,total generation rate is a 
monotonically decreasing function of depth. In other words, the peak 
carrier generation for a homojunction cell occurs. at the front surface. 

However, the carrier generation distribution for a 1 variable composi- 
tion cell has a more controllable shape determined by the composition 
profile and the wavelength content of the .solar spectrum , As can be 
seen from Figure 2.8, grading from a semiconductor with ,a large bandgap 
to one with a smaller gap. can lead to a generation rate peak well below 
the surface of the cell. As the composition is varied so as to reduce 
the bandgap, a larger share of the incident photons have sufficient 
energy to generate carriers. Figure 2.9 shows the extreme case of an 
abrupt change from a wide bandgap material (AlAs) to a much smaller one 
(GaAs),. The AlAs is relatively transparent to most of the solar photons, 
while a large percentage of these photons are absorbed by the GaAs layer 
beneath. This "window effect' 1 is a major advantage of variable 
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Figure 2.8 Generation rate vs. depth below surface for a graded A1 Ga _ 1 
solar cell (Figure 2,6&) under AMO illumination with an x -x 
0.07uM antireflection layer. 
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composition solar cells because it separaxes the region of maximum 
carrier generation from the top layer of high surface recombination 
velocity associated with surface imperfections. In addition, the 
transparency of the window means that the surface layer can be 
relatively thick and thereby reduce series resistance losses that 
plague GaAs solar cells [58]. Of course, such a configuration also 
introduces interface states that reduce lifetime at the heterojunction 
so that a poor lattice match between materials can nullify the gain in - 
collection efficiency that would otherwise be brought about by the 
window effect. The problem of interface states is discussed further 
in the next section. 

2.2.3 Interface Recombination 

Although it is to be expected that carrier loss by surface 
recombination can be effectively controlled by use of a wide bandgap 
window layer, -it is by no means certain that interface states due to 
dangling bonds at the hetero junction will not prevent an overall 
improvement in collection efficiency. Energy conversion efficiency is 
especially sensitive to interface state density when the heterojunction 
is placed within the- p-n junction depletion region because recombination 
there can significantly increase the diode dark current. If depletion 
region recombination is large enough so that electron and hole currents 
cannot be assumed constant across the depletion layer , then the current 
density may be more nearly proportional to exp(qV/2kT) than to 
exp(qV/kT) given by the ideal diode law. In other words, interface - 
states within the depletion region can soften the ”knee u of the J-V 
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characteristics and thus 1 reduce peak output power by reducing the .curve 
factor as well as the open circuit voltage. 

Howeyer, ther.e are several design options available to, control 
losses due to interface recombination. The most obvious is the .selec- 
tion of a pair, of semiconductor materials with a close lattice match., 
For example, an AlAs-GaAs hetero junction, with a lattice mismatch of 
about 0.14%, would appear to be a likely candidate for minimizing 
interface states. An interface between GaP and GaAs, on the other hand, 
gives a lattice mismatch of about 4%, and could be expected to suffer 
far more from interface recombination losses. 

A second method of minimizing the effects of interface recombina- 
tion is to place the heterojunction away from the center ,of the p-n. 
junction depletion region-. This prevents the drastic increase in dark - 
current previously mentioned, and also removes the recombination sites • 
from the area of maximum minority carrier collection due to the built- 
in field within the depletion region, ©f course, some loss will still 
occur- as carriers recombine at the ^dangling bonds even when they are 
placed in the quasi-neutral region of the device. In order to further 
reduce the effect of recombination at lattice mismatch sites,- the change - 
from the wide bandgap material to the smaller bandgap material can.be 
spread over a finite distance until the effective lifetime in the bulk 
is- not significantly less than the bulk lifetime without interface 
states . 

The detailed .reasoning behind these ideas concerning interface 
recombination is discussed further in Chapter 4 which deals “with the. 
computer, modeling of several device phenomena. In addition, Chapter ,6 
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discusses the results of several computer simulation runs made expressly 
to evaluate the effect of interface states on various material arrange- 
ments . 

2.2.4 Optimization of the Bandgap in the Substrate 

It was mentioned in Section 2.1 that no single material can have 
a bandgap such that all solar photons have sufficient energy to excite 
carriers without wasting any excess energy as heat. Aside from the 
secondary generation of pairs by some highly excited carriers , any 
energy relinquished by absorbed photons exceeding that necessary to 
bridge the bandgap, eventually serves only to heat up the lattice. 

On the other hand, virtually all those photons with energy less than 
the gap are lost as they pass through the crystal without being absorbed 
in significant numbers. 

Research by others indicates that the optimum compromise between 
these two loss mechanisms calls for a material with a bandgap of about 
*1.4eV [4,21,23,54]. Although GaAs has a bandgap of approximately this 
value, it cannot be expected that a single semiconductor will happen 
to have the optimum gap as well as the most desirable electrical 
characteristics. It would be more reasonable to examine an alloy 
between two semiconductors with the proportions selected so as to give 
optimum overall conversion efficiency. The results of computer simula- 
tion runs to determine the optimum composition beneath the window are 
presented in Chapter 6. 

Another option to improve absorption efficiency is also available 
to the device designer. If composition is graded so that light passes 
through layers of progressively smaller bandgap, photons of 
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progressively lower energy can be picked off with less energy waste 
than would be possible with any fixed bandgap . However, given the com- 
plex profile of the solar spectrum, and the fact that the bandgap 
cannot -be controlled independently of other material parameters, it is 
clear that the determination of -the composition profile for maximum 
conversion efficiency is not a simple problem . In such a situation, 
the computer simulation daxa, to be discussed in Chapters 5 and 6, can 
be very useful. 


2.2.5 Built-in Fields due to -Composition Grading 

As mentioned in Section 2.1, homogenous solar 'cells depend on 
doping gradients for the creation of built-in elecxric fields that 
collect the carriers generated by optical absorption. However, when 
material composition varies with position, new built-in electric fields 
are established that can be used to improve carrier collection 
efficiency . 

The effective built-in field that operates on electrons is 
proportional to the slope of the conduction band edge as 


Similarly for holes, 
valence band edge as 


no q 


dx 


(2.7) 


the effective field depends en the slope of the 


E 

po 


, dCEp-E, ) 

1 F y 

q_ dx 


( 2 . 8 ) 


since 
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(2.9) 
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It is clear that different effective fields can he made to operate on 
electrons and holes only if the bandgap E is a function of position. 
The difference in the two fields is 


dE 


E — E = - 
no po ~ q dx 


( 2 . 10 ) 


Therefore, by varying the bandgap with position, a built-in field can 
be established to help sweep excess minority carriers to the p-n 


junction and thus improve collection efficiency. This option is another 
advantage gained by allowing material composition to vary with position. 
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3 . DEVICE EQUATIONS FOR COMPUTER ANALYSIS 
3.1 Introduction 

In order to adequately analyze variable composition sqlar cells, 
several device phenomena must be incorporated into the mathematical 
equations used to model the cells. Among the effects important in 
ordinary homojunction solar cells are: 

» 

l t Drift and diffusion currents 

2. Position dependent doping 

3. Doping dependent mobility 

4. Optical carrier generation 

5. Bulk generation-recombination effects 

6. Surface recombination effects 

The introduction of spatially varying composition demands that 
additional factors be accounted for: 

1. Position dependent bandgap 

2. Position dependent electron affinity 

3. Built-in fields due to a varying bandgap 

4. Composition dependent refractive index at the surface 

5. He ter o junction interface recombination 

6. Other position dependent material parameters such as: 

a . mobility 

b. dielectric constant 

c. optical absorption coefficient 

The computer program used to simulate solar cell behaviqr is the 
result of work by several individuals over a number of years. In its 
original form it was used to analyze p-n junction diodes with no 
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provision for .optical carrier generation -[13] . In 1974 it was .modified 
by introduction of -a separately computed optical generation rate so that 
homogeneous- solar cells could be studied,. C 19 3 • In r 1975, as part of the, 

present work, the program was again modified to allow for position 
dependent composition * In the following sections, the* equations used 
in the homogeneous solar cell program .will be presented, followed -by a 
development of the corresponding equations used to model variable ' 
composition solar cell .behavior 


3.2 Equations for Homogeneous Solar Cells 
The one-dimensional .device equations used to describe homo junction 
solar cell operation are 
Transport: 


Continuity: 


, _ _ 9n 

J = qu nE + q D — , 
n * n u n 9x 


(3.1) 


J = qjj pE 
P P 

u p 9x 

(3.2) 

|a,. u + 

i 3J « 

± E + n 

q 9x e’ 

(3.3) 


( 



, 3J 

-r-E-t G , 
q 9x e 

(3.4) 


Poisson’s equation: 


||= f (p-ntN)i 


Auxilliary equations: 


np-n. 


U = 


t (p+p, )+t (n+n n ) : 

no r *1 po 1 


(3.5) 


(3.6) 
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( Shockley-Read-Hall model with single trapping level) 

Ep -E 

n = N c exp[" -^y- C | 3 (non-rdegensrate) , 

E -Ep 

p = exp[ — ~ (non-degenerate)-,. 

• N = - N a (ne.t -doping) 


ta.7) - 

(3.8) 

(3.9) - 


This system of equations is valid only -for a 'device made of a 
single semiconductor material doped with enough impurities to- create 
the desired doping profile.- The electron energy- hand structure and its 
associated characteristics must-be the same throughout the device. 

As shown. by Graham and Dunbar [13 ,19] , the .preceding group of 
equations can.be reduced to three equations in the three unknowns (j> n , 
and <j>p, where $ is the electrostatic potential and <j> n and $ are the 
Fermi -potentials defined -by 

<j) = - - E - , (3.10) 

n q Fn- 

<|> =, - i E_ . (3.11) 

p q Fp 


It can be shown thaf the same basic approach using the same ’three 
variables,, can.be used to, analyze variable composition solar cells. The . 
modifications required to account 'fot material variations are discussed 
in the -next section. 


3 » 3 Modification of Equations to Allow 
for Variable ■Composition 

In order to rewrite the device equations so that they apply .to 
variable .composition devices, it is useful to recall the basic principles 
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from which .they were derived . In particular » Equations ' (3 .1) and (3.2) 
can be written in the more general form 


J = 
n 


dE 


U n 
n 


Fn 


dx 


(3.12) 


J 

P 


dE. 


= V 




dx 


(3.13) 


Transport • equations involving the field ’and, carrier concentration 
gradients proceed from the above equations as follows. . -Referring to 
Figure 3.1 the conduction band edge energy can be. written in, general as 


E c = “ V (3.14) 

where x is the electron affinity, if) is the electrostatic potential and 
c 

ifi Q is soma 'ref erence potential. In a similar manner 

E = q(i|> -tfO-x -E (3.15)- 

V ^ r o eg 

From Equations (3.7) and (3.8), 

E,, = E + kT in (--), (3.16) 

Fn c N 


E = E v - kT £n (E-). (3.17) 

P V 

Then using Equations (3.14) and (3.16.), Equation (3.12) can.be rewritten 
as 


dx 

J s p n CqE - , 
n n ^ dx 


dN 

c + Ml _ Ml Q- i 

• n dx ~ K "dx 
c ■ 


(3.18) 


In a similar manner 'the. hole density equation becomes 

5k _ 5s. .. + kT - dNy i 

dx dx p dx N dx J 
^ v 


J p = U p p [qE - 


(3.19) 
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Figure 3 0 1 Elec-cron energy band diagram for a heterojunction solar 
cell at equilibrium in the dark* 
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For a homogeneous ' material , x q -s E g s N^, ;and are .constant through- 
out the device, and Equations (3.18) and (3.19} simplify to Equations 
( 3 . 1 ) and ( 3 „ 2 ) . But , for a , variable ■ composition structure , Equations 
(3.1) and (3.2) musx be replaced by Equations , (3.18) and, (3.19), or by. 
the simpler forms of Equations (3.12) and (3.13). As discussed in 
Chapter 2, effective fields acting on electrons and holes .arise in 
Equations (3.18) and (3.19) from bandgap, electron affinity, .and, density 
of .states variations with position., 

The (continuity Equations (3.3) and (3.40 apply ^equally well to 
heterogeneous or homogeneous materials i Poisson *s .equation, however,, 
fails to .reduce 'to t the simple expression of Equation (3.5) when composi- 
tion (and ’therefore dielectric constant, e) is allowed to vary with 
position. The more general ‘foqm whiph results from v V • D = p, is 



S* (p-n+N) - 


1 dip de 
e dx dx° 


(3.20) 


Thus,- when the dielectric .constant becomes a .function ,©f position. 
Equation (3.20) must replace Equation (3.5). 

The -carrier concentration equations demand.no modification,- but 
they can be written in a more convenient form for the purposes of com- 
puter-analysis. Using Equations (3,10), (3,11), (3.14)- and (3,15),- 
Equations ' (3.7) and (3,8) become: 


n = ,n i2 expCja (*-* n t — - * c + — tn.j-j)] 


■kT 


(3.21) 


X E vt N 

p = n, exp Cr§- (<j> -i|) ^ + ip t — An -^*-)J 

r x2 ^ kT p q q o q n^ 2 


(3.22) 
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where n.~ = intrinsic carrier concentration of some reference material 

i2 

(taken as alloy component: #2 in this work). 

The composition dependent (and therefore position dependent) para- 
meters, can , be conveniently grouped-, into two terms 


6 ■ 
n 




(3.23) 


N 

e = - - (JC +E_) + .r + —in (— -) 

P q c g T o q n i2 


(3.24) 


Then, the carrier density equations can be written concisely as 


n = n. 2 exp (<M n + 


(3.25) 


✓ 

p - n i2 exp C 4>p“ ) 3 - (3.26) 

These expressions are analogous to the forms used, in the computer 
analysis of homogeneous solar cells [193, i.e. 

» = n i exp c i# 'f-y 3 - 

P = n. exp (<f» -tj>)3. 

The parameters, 9^ and 9^ depend only on composition, and. must be 
externally supplied as functions of position for the specific composition 
profile being considered. 

The other auxiiliary equations remain the same whether or not com- 
position is spatially varying.. However, it is important to note that 
several variables no longer can.be assumed to be constant. In fact, the, 
material .parameters u , y , e, n., n. and p. must be given as known ,func- 

H p 1 * X «I» 

tions of composition (or .position, when the -composition profile is. known). 
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Thus, the> complete one-dimensional mathematical ’model for -variable 
composition solar cells consists of The following system of equations: 


d<f> 

J =-qu n ~r — , 
n ^ n dx 

(3.27) 

d<j> 

J = -qu p 
p ^ p dx 

(3.28) 

. , dJ 

9n _ p n + 1 ‘ n 

a.t " G e ° + q dx 9 

(3.29) 
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3t e U q dx ’ 

(3.30) 

ik . _ a. 

,2 e e dx dx 

dx 

(3.31) 

np-n. 2 

*u = — — 

t IJ0 <ptp i ) + ' T po (p+n i ) . 

(3.32) 

n = n. 2 exp (*-* n +8 n ):. 

(3.33) 

P = n 12 exp 

(3.34) 


with the, following parameters given as functions of position: 

N - N(x) = - N^, net doping profile (3.35) 

C = C(x),' composition profile (3.36) 

(for example, G.= mole fraction, x, of 

AlAs In an A1 Ga.. As solar cell) 
x 1-x 

V - y (x), electron mobility (3.37) 

n n 

Up - Up(x), hole mobility (3.38) 

e = e(x), dielectric constant 


(3.39) 



0& *« a 


n i ~ ’ intrinsic carrier concentration 


(3.40) 


n^(x), trapping center parameter for electrons (3,41) 


p^(x) , trapping center parameter for holes 


( 3,o 42 ) 


T (x) s excess carrier lifetime in p+ material (3.43) 
no 

x (x) 9 excess carrier lifetime in n material (3,-44) 
po 


E = E (x) , bandgap 
g g 

6^ = Q (x), conduction band parameter 


(3.45) 


(3.46) 


9 = 9 (x), valence band parameter 

P P 

G = G (x), optical generation rate 
6 0 


(3.47) 


(3,48) 


The position dependence of these parameters can arise from compositional- 
changes and/or doping changes with position. The parameters e, 8^, 0^, _ 
n^, n^j are fundamentally functions of composition .C(x)„' The param- 
eters and are functions or both doping density and alloy composi- 
tion. The generation rate, G , has an intrinsic position dependence as 
well as a composition dependence. Finally and are assumed to 
depend on the density of intrinsic recombination sites as well as the 
density of dislocations introduced, by lattice mismatch in graded, compo- 
sitional materials. Magnetic fields and thermal gradients are assumed 
to be negligible. Additional discussion of the modeling of these various 
material parameters is given in Chapter 4. 

This apparently unwieldy mathematical system can actually be reduced 


to three, equations in three unknowns in exactly the .same manner used for 
homogeneous solar cells in previous computer models. The fact that 
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composition (and the associated material parameters ) is allowed to vary 
through the, cell simply adds some complexity to the detailed computations 
without altering the method of solution* The next section discusses the 
techniques used to solve the device equations by computer. 


3.4 Computer Solution by Quasilinearization 
Although previous reports by Graham and Dunbar £13,19] will prove 
useful as references for the detailed computer techniques ,- this section 
discusses those ideas concerning the numerical solution of the .device 
equations that are most important in understanding the. analysis results 
that follow. In order to maximize the accuracy and minimize routine 
calculations, all basic variables are normalized by the constants given 
in .Table 3.-1 before the solution process begins. The normalized steady 
state device equations for variable composition solar cells are, from 
Section 3,3 and Table 3.1, 
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Table 3.1 Normalization constants (Subscript "2 11 indicates semiconduc- 
tor #2 in alloy, GaAs in this work) . 



Normalized 

Variable 

Description 

Normalizing 

Constant 

X 

position 

* 

L D = /e 2 V T /qn i2 

T ' ,T 
no po 

lifetimes 

L D /D o 

<P 

electrostatic potential 

V T = kT/q 

<j> , <j> 
n p 

quasi-Fermi potentials 

V T 

V A 

terminal voltage 

V T 

E 

electric field 

v t /l d 

n »P 9 n jyPj_ 

carrier densities 

n i2 = /!i c2 N v2 e!tp( - E g 2 /ki,) 

n,n d ,n a 

net impurity, donor, 
acceptor concentrations 

n i2 

J ’ J n> J p 

current densities 

^ D o n i2 /I 'D • 

D ,D 
n* P 

carrier diffusion constants 

2 

D = 1cm /sec 
o 

j V 

n p 

carrier mobilities 

V V T 

U 

generation-recombination rate 

D o n i2 /L D 

G 

e 

optical generation rate 

Vi2 /L D 

s 

recombination velocity 

D o /L D 

e ,0 
n’ p 

band structure parameters 

V T 

N ,N 
c’ V 

E 

g 

effective densities of states 
bandgap 

n i2 , r «-a 

X c 

electron affinity 

kT 

£ 

relative dielectric constant 

c 2 
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r no (p+P l ) * T po (n+n l ) 


(3.56) 


Taking if) , <f> 5 and cf> as The sex ef independent variables , xhe above 
n P 

eight equations can be reduced to three differential equations ef the 
form: 


TT* F 1 C *> V r) 

dx r 


(3.57) 


d 2 <f> 


n 


dx' 


2 =f 2 <*, * n9 ^ p5 *•, ♦;) 


(3.58) 


d 2 <p 

F 3 ( V V *'* 
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where the primes indicate differentiation with respect to x. For refer- 
ence, xhe three functions F^, F^, and F^ are given below. 

F l^’ *p # ^ = e* £exp(^)-<f) n te n ) - exp(<j> -if)+0 ) - N - ~ ||-3 (3.60) 
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As discussed in Chapter* 4, the parameters e, G^, N, n^, p^, T no , t^ q , 

0 . 0 » u » U must be evaluated as functions of position, 
n/ p ? n\ p 

The method of solution used in the present research is identical to 

that used by Graham and Dupbar for less complicated structures [13,19]. 

Briefly, the computer program accepts input data specifying device 

structure and terminal voltage, and uses first order approximations 

and/or available data on previous ' solutions for a similar device to 

produce an initial estimate for each of the variables in the arguments 

of F^, and Fg. Linear .approximations to Equations (3.57-” 3,59) are 

then obtained by expanding F^, and in function space about the 

initial estimate and truncating the resulting series to eliminate higher 

order (nonlinear) terms, These, approximations are used to generate a 

new estimate for the variables ij>, <j> n and $ and the process is repeated.,. 

If the initial estimate was close enough to the exact solution, this 

iterative procedure will converge to a solution of the original set of 

equations iTe . a position dependence for <j> , and 4> throughout the 

n p 

solar cell. Further calculations on these three variables provide pro- 
files of physically meaningful quantities' such as electron and hole . 
concentrations and current densities everywhere within the device. These 
results can then be examined, to insure that the program's iterative 
process has converged to a physically reasonable solution. 

The operation of ’the main solar cell analysis program is summarized 
in Figure 3.2. First * the material parameters for the ternary alloy 
system are initialized and the specifications (doping levels, layer 
thicknesses, polarity, etc.) for the desired structure are read in. 

Th,e terminal voltage is then read in, and the program calculates an 




Figure 3.2 Flow chart for main variable composition solar cell analysis 
program. 
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inixial estimate for the 1 three independent variables ip, <J> n# and . If 
the same device structure has been previously analyzed (perhaps at 
slightly different lighting and voltage levels) data may be recalled from 
disk storage to provide an initial estimate. However, if no previous 
analysis has been done, .the program must 'build up ap initial estimate 
in stages, beginning with an analysis of a homogeneous , cell in darkness 
with zero volts across the .terminals . The - compos ition profile is then 
built up in 'Steps by using the previous solution as an initial estimate 
for dexermining a more exact solution by 'quasilinearization. Once the . 
composition profile is established, the terminal voltage is stepped up 
in a similar fashion until the 'expected open 1 circuit voltage has been 
exceeded. Next, the illumination intensity (i’.e., generation rate) is , 
built up in steps. As soon as this stepping process is complete, the 
terminal voltage can be, swept down again 'to obtain a complete J-V 
characteristic for the cell nrjder full lighted ,cpnditions. 

As shown in the flow chart,, the .results can 'be obtained .in both 
plotted and tabulated form. Figures 3.3 through 3*12 present a sample 
of the plotted output xhat the analysis program produced_ for an AlAs-> 

GaAs he.tero junction solar cell with a window thickness of 0.1 yM and 
an abrupt -p-n junction placed 0.4 yM below the surface. The abrupt 
change in bandgap and, electron affinity, at the '0.1 pM he tero junction 
is clearly evident in the ‘equilibrium band diagram^ Figure 3.3. The 
built-in fields 'around the hetero junction and p-^n ^junction can be seen 
in Figure 3.4, which shows the .electrostatic potential near the surface 
under AM© -conditions at various terminal voltages.. Figures 3.5 and 3,7 
show the position dependence, of electron and hole .current densities for 
the entire solar cell depth while Figures ,3,6 and 3.8 are, expanded .views 
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Figure 3 0 5 Electron current density through the entire depth of an 
abrupt heterojunction solar cell. 
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Expanded plot of electron current density near the front 
surface , 
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Figure 3.12 Dependence of efficiency on terminal voltage for an 
abrupt heterojunction AlAs-on-GaAs solar cell as 
calculated by computer simulation. 
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of the same data near the front surface. Although the carrier densities 
plotted in Figure 3.9 appear relatively smooth, an expanded view 
(Figure 3.10) near the surface shows the large changes in carrier 
concentration that occur around an abrupt AlAs-GaAs hetero junction and 
an abrupt p-n junction. 

The J-V characteristic calculated by the computer analysis program 
for this device is shown in -Figure 3.11, and finally the dependence of 
peak efficiency on terminal voltage is plotted in Figure 3.12. The 
availability of data such as that shown in the preceding ten figures 
makes computer simulation a useful tool in determining the device 
characteristics that control solar cell behavior. 
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4. MATERIAL PARAMETER MODELING 
4.1 Introduction 

As mentioned in earlier chapters, various material parameters must' 
be given as functions of position before the .device equations can be 
solved. The material parameter modeling discussed in the following 
sections makes possible the computer analysis of solar cells made of 
a spatially varying alloy of two compatible semiconductors. Given the 
material properties of the two binary materials , the_ computer program 
generates approximate parameters for the specified ternary composition 
profile. For example, solar cells made of the alloy, A1 Ga- As, are 
analyzed, by providing the material properties of AlAs and GaAs , along 
with the desired composition profile, ,i:e_« the mole fraction of AlAs, x, 
vs„ position. The program then uses this information to determine the 
values of the parameters needed for ,the device equations at each point 
within the solar cell. Equations of both a .theoretical and empirical 
nature have been used to obtain the alloy material parameters from the 
corresponding binary material parameters. . Whenever possible, these 
approximations have been checked against experimental results and - 
adjusted to improve the .agreement with available data. The most impor- 
tant approximations ' are described- in .the following sections. 

4.2 Dielectric Constant' vs. Composition 
Given the dielectric constants, e 1 ,and e 2 , for semiconductors 1 and 
2, the following interpolation scheme is used to estimate the dielectric 
constant, e, for. an alloy that has mole fraction, C, of material 1, 

[ 24 ]: 
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e 


El -1 e.-i 

l+2tC(~r) + (l-C)(-~r)J 


‘V 2 


V 2 


s -l 
1 x 


- 1 C( ’e 1 +2 :) “ ( ' 1 C ^e 2 +2^ 


(4.1) 


This form was used to approximate the values of the high frequency- 
dielectric constant, e^ 9 and the low frequency dielectric constant 
£j> , both of which are needed in the mobility calculations to be 
discussed in, later sections. 


4 o 3 Band Structure Parameters vs. Compose trion 
4„3„1 Introduction 

A1 Ga. As and GaAs P have both a direct and an indirect bandgap 

X J."*X J.”X x 

which are important over various alloy compositions, and the solar cell 

program permits the specification of two valleys, as shown in Figure 4.1, 

each with its own composition dependent parameters such as bandgap, 

effective mass and mobility. Since the mean time between intervalley 

—12 

scattering is about 10 sec (or less), these two electron populations 
are essentially at equilibrium with each other and a single electron 
quasi-Fermi potential, <j>: , can be used in the device equations. Then an 

II i 

"effective" electron mobility and "effective" electron mass can be 
defined such that a single current density equation (Equation 3,27) and 
a single carrier density equation (Equation 3,33) is sufficient for this 
two conduction band valley model. The following sections describe the 
techniques used to determine the resultant band -structure parameters 
(such as electron and hole mobilities) of an alloy, from. the band 
structure parameters for each of the two components „ 

ORIGINAL WG* 

OB! POOR QUAU3K 




Figure 4oi Electron energy. E, vs» wave vector , k, for a) indirect 
semiconductor, material 1, b) direct semiconductor, 
material 2, c) alloy of materials 1 and 2, 
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4o3„2 Bandgap^ vs . Composition 

Experimental studies have shown that the direct energy gap for 
III-V alloys has an approximately quadratic dependence on composition 
such as 


E , = aC 2 + bC + E (4.2) 

gd gd2 

where a and b are experimental values,, The same form was used to model 

the indirect variation in GaAs., P , while for A3. Ga. As, the indirect 

1-x X s x 1-x 

gap was estimated by a linear variation with composition. The bandgap 
is the minimum of the direct or indirect bandgaps . 


4„3<,3 Effective Masses 

The next task is to determine the effective masses for holes and 
electrons. Since the valence band in each of the two semiconductors is 
assumed to have a single dominant valley with a characteristic effective 
mass, it is simply necessary to interpolate between the two effective 
masses to approximate the resultant hole effective mass for the alley. 
The following interpolation form was used [24J: 


m 


m- 


Pi 


T — X 

irr 


1-C 


(4c 3) 


P 2 


where m* , = hole' effective mass in material 1 
pi 

m* _ = hole effective mass in material 2 

p2 

m* = hole effective mass in all oy 
P J 

C = mole fraction of material 1 

The same form was used to find the electron effective mass for rhe 
indirect valley: 



C 1-C 

m* , _ m* 

cdl cd2 


(4.4) 
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where conduction band effective mass in direct valley for 

material 1. 

m * c £2 = conduction band effective mass in direct -valley for 
material 2. 

m* cd = conduction band effective mass in direct valley for alloy. 
Likewise for the indirect valley: 


C , 1-C 

- + -x 


iri" . m- . , m 5 . _ 
ci cil ci2 


(4,5) 


Now, it will prove useful to define, if possible, a single effec- 
tive mass for all electrons. This can be done by considering the total 
electron .density expressed as 


» * "=d + N ci 


(4.6) 


where N ^ = direct valley effective density of states t 
N ^ - indirect valley effective density of states 
E c £ = direct valley minimum 
E . - indirect valley minimum 


It is desired to obtain the following form; ‘ 


n “ 


(4.7) 


where = effective density of states for entire conduction band 

E - absolute minimum of conduction band (the' lesser of E , and 
c cd 

E 0i>- 

Si „- N cd * K “"'cd 3 ' 2 

u _ v a 3/ 2 
Ci ci 

N - K m* 3/2 
c n 
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where K is a proportionality constant and - effective mass for all 
electrons, then 


, 3/2 r E Fn~ E c n * 3/2 r E Fn E cd-, * 3/2 r E Fn" E ci 

' expC— pfr — 1 = exp[— ^ ] + m* . exp[— 


n 


m % = «od 3/2 ex P ( -% 2i) + 


kT 

3/2 


■1 


/ E c E ci^2/3 
mw ci expC-^r-)] 


or, in terms of bandgaps: 


m *n = Cm<, od 3/2 * ^i 372 ex P (-£ifSi)3 2/3 


(4.8) 


where E , , E ,, E , m* . and m* . are given by Equations (4.2), (4.4), 
gd* gi 5 g * cd ci 

and (4.5). 


4.3.4 Hole Mobility 

Mobility depends not only on alloy composition, but on the doping 
lqvpl as well. In order to calculate hole mobility for an arbitrary 
doping and composition, an empirical formula, V was constructed to 
apprpximate the doping dependence of hole mobility for one of the alloy 
components , to be referred to as material 2: 


v = V 10 


wherq N = total impurity concentration. The form of f^ was taken to 

be 



(N) 


A 

£ . 

(1+N/N )9 
P 


(4.9) 
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Npw, because 


q<x > m 5; „ < x > 

u = P. = u _Pi_ _P_ 

m* p2 <t _> m* s 

P p2 p 


(4.10) 


y/here <x > and <x > are the mean times between scattering for holes in 

f P 

the alloy apd in material 2, respectively, it is apparent that if 

ft * and pan bp specified as functions of composition, then a 

P “ P 

p^asppable approximation to hole mobility for any doping level and 
coiJtpogijtiop cap bp made. In the III-V semiconductors, it was assumed 
{nobility, ip controlled primarily by a polar optical phonon 
spattering process 1 of the form £253 

K 


<T> 


/ m*" (- 


1 


•) 


(4.11) 


K ip 'p proportionality constant for all materials and and 
ppe tljp fiigh and low frequency relative dielectric constants. 

£qmbiping Equations '(4.10) and (4.11) then gives as an approximation 

tp mobility 


Mp(N»P) = 


f p2 (N) 


i* 3//2 ( 

P 2 ^ 2 






(4.12) 


“V 

V y ? R d ro<i p 


^nd e h2 are known parameters of material 2, and 
are given as functions of composition by. Equations (4.1) 


3r$ (4,3). 


7TT 


y aridity of this assumption has been questioned for hole mobility 
in p-type indirect gap and p~type III-V semiconductors [48], 
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4^,5, Electron Mobility 

The modeling' of electron mobility is more complex than hole mobility 
because pf the indirect 'and direct bands. Treating the direct and 
indirect valleys independently, and using the same Technique as described 
£o.r t]ole mobility, a direct valley electron mobility, y^, and an indirect 
valley electron mobility, y. 9 can be determined approximately as: 


y d (N,C) = 


y i (N',C) = 


f (N) -m* 3/2 (~- - •—-) 

n2 Cd2 e h9 

3/2 f 1 lv 
od E h % 

f (N)m* 3 ^ 2 (— _L_) 

nl Gi 2 ---- E -h-l Sui- 

3/2, 1 1. 

m- . (- — J 


(4.13) 


(4.14) 


where m$ cd2 » m *ail* \ 2 5 e £l s s h2 and e hl aX ' e known P ara ^eters of 
{Hatqpialp 1 and 2, f (N) is an empirical function relating electron 
mobility to doping for material 2, f^(N) is an empirical function 
plating electron mobility to doping for material 1 and m* cd and 

ppe given by Equations (4.1), (4.4) and (4.5). The functionsj.f^ and 
f us^d the same functional form as Equation (4.9). An "effective" 
mphi-li^y for use in the current density equation can be found by 
weighing th® direct and indirect mobilities by their respective electron 


pppqlatiops as 


U n .('%C>. ®.yr d R£ + y i (l-R d ). 


(4.15) 


whppp is the fraction of electrons in the direct valley given by the 


expression 


R = ■ ■ 
d n,+n. 


,,+n. m* . E ,-E 

d 1 ' i + (- r gA) 3 / 2 'exp(- S-t-SA) 

cd 


(4.16) 
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,4,3.6 Electron Affinity 

It is apparent from experimental data reported by others [263 that 
the electron affinity for a direct gap ternary alloy of two binary 
semiconductors is often approximately linearly related to the bandgap as 


X c " x c2 " K a ^ E gd " E g2 ) 


(4.17) 


where x c 2 the electron affinity for one binary semiconductor (mater- 
ial #2, having a direct bandgap), E is the bandgap of material #2, 
is the direct bandgap of the alloy, and is a proportionality 
ppnstant associated with the alloy being considered. 

The determination of electron affinity is slightly more complicated 
if one binary semiconductor (material #1) is indirect so that the alloy 
itsplf will be indirect over some range of composition. Since electron 
affinity is defined as the energy difference between the conduction band 
edge and the vacuum level, the relationship between absolute bandgap 
and electron affinity in the present work was approximated by 


X = Xo~K(Ej-E 0 )+E.-E 
c . A c2 a gd g2 gd g 


(4.18) 


wfiere E is the absolute gap of the alloy. This equation preserves the 
linpar relationship between the 'vacuum level and the direct valley 
conduction band minimum over both direct and indirect alloy compositions. 


4,3.7 Band Parameters, 0 and 0 
> n 

In order to model 'the composition dependence of the band parameters, 
0^ and 0^, defined in Section 3.3, some knowledge of the variation of 
electron affinity, x a and bandgap, -£ , must be available. It has been 
fopnd to be advantageous to select the potential reference, so that 
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6=6=0 when the mole fraction of material 1 is -zero . In other 
P P 

words , 0 n and 9^ are zero in regions of- the solar cell consisting 
entirely of material 2. Then, from Equations (3.23) and (3.24) 


. a c2 kT . / c2. 

^ q q n. 2 .’ 


(4,19). 


pr equivalently, 


X„„ +E 


«!> = 
o 


S2 _« tn ( .V2>, 


(4-. 20) 


NqWs inserting Equation (4.19) into Equation (3.23), and Equation (4.20) 
iqto Equation (3.24) gives 


X C “X C 2 v-j G 

0 = _£_££ + H. «, n ( -2^.) , 
n <1 <1 N 0 


(4.21) 




fi^lly, in terms of effective mass 


„ c , O Si . / IK 

9 = + — — in (-5T-) 

n q , 2 q m* 2 


-(Ay +AE ) 0 . _ m* 

t 3*l, n (_E.) 


w]ier„e AE^ = E - E g2 s ~ X c ~ X c2 are gi ven as functions of 

qomposition by Equations (4.2) and (4.13). - 

Equations (4.23) and (4.24) are the forms used to compute the band 
parameter^ for ,an arbitrary ternary alloy of two .binary materials . 
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4,4 Interface Recombination 


At an abrupt interface between two semiconductors, lattice mismatch 
introduces trapping levels due to dangling bonds . The density of 


cjapgling bonds at such an abrupt hetero junction is given by [273: 

2 2 , 


N h = 


K I I 
2 2 
a 2 a i 


(4.25) 


where and Are the lattice constants of the two materials and K = 4 
for a [1003 interface, K = 2/2 for a [1103 interface, and K = 4// 3 for 
a £lll>3 interface. If the composition change is gradual, rather than 
abrupt, the trapping centers have been assumed to have a volume density 
that is proportional to the slope of the composition profile as 


N = N ~ 

g h dx 


(4.26) 


These interface states have been incorporated into the Shockley-Read-Hall 
model by modification of the excess carrier lifetimes: 
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x 

no 


no 1 


Bulk 


+ I c 
g R 


(4.27) 
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T 

Pe 


1 

p0 I Bulk 


NO, 

g R 


(4.28) 


where is the capture coefficient. Etxehberg and Kressel £283 obtained 
3 

a value of 8x10 cm/s for the interface recombination velocity at an 
fifopupt A1 q 5 GAq gAs - GaAs hetero junction, from which C R can be calcula- 
ted as 
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0.5M, 


h _ 


4 l a 2 ‘ a ll 

(Q.5)[— 

a 2 a l 

3 9 

8x10 cm/s 


C ='5,2 x 10 8 cm 3 /s, 
K 


where the lattice constants and for AlAs and GaAs were taken as 
5.661 x 10 8 cm and 5.654 x 10 8 cm, respectively [293. 


4.5 Optical Carrier Generation 

It is apparent from the equations discussed in Chapter 3 that an 
optical generation rate G g for electron and hole pairs must be known as 
a function of position prior to the solution of the equations. This 
generation rate profile is calculated by means of a separate computer 
program in a manner similar to that used for homogeneous solar cells 
[133. 

The program calculates the position dependence of the generation 
rate for a device made of an alloy of two compatible semiconductors 
covered by an antireflection layer of arbitrary thickness and made of 
SiO (see Figure 4.2). 

The composition profile is arbitrary, and the program is designed 
to determine the most efficient thickness for the antireflection layer. 
The input data required consists of the optical properties (absorption 
coefficients and refractive indices) of the two semiconductors and SiO, 
and the intensity vs. wavelength profile of the solar spectrum. 

Although only AMO conditions have bpen selected for calculations 
reported in this work, data for AM2 conditions is also available as a 
user option. The intensity profile of the solar spectrum is provided 
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SiO 

Antiref lec cion 
Layer 


Figure b -.2 Device sxrucxure assumed for calculaxion of generation 
rate profile o 
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in terms of power density per wavelength interval P(A), at a finite 
number of wavelength values. This information then gives the incident 
photon density per second as 

N. (A) = - P p - >vvv (4.29) 

l E(A ) 

where E(A) is the average photon energy over the wavelength interval AA. 

The program uses the absorption coefficient and refractive index 
data for SiO and the semiconductor alloy to determine a transmission 
coefficient T(A) for the antireflection layer. Thus the photon density 
per second entering the semiconductor (at ,x = 0) is * 

N(AjO') = N.U) T(A) (4.30) 

At this point however, the calculations for variable composition 
solar cells begin to differ from those for homogeneous cells because 
the absorption coefficient a is a function of position as well as 
wavelength. If N(A,x) is the photon density per second passing any 
point at depth x beneath the semiconductor crystal surface, then the 
rate of change of photon density can be written 

~^=-.a[A,C(x)J NU,x) (4.31) 

where C(x) is the “mole fraction of binary semiconductor #1 (such as 
AlAs) present in the ternary alloy (A1 Ga. As for example). 

Since experimental absorption data is available for only a limited 
number of alloy compositions, an interpolation procedure was created 
to determine a for arbitrary wavelength and composition. Because the 
program allows for a two-conduction-band-valley model for an alloy 
between an indirecx gap material (#1) and a direct gap material (#2) the 
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total absorption coefficient for a given wavelength was taken to be the 

sum of an absorption coefficient due to direct transitions and an 

absorption coefficient due to indirect transitions. The indirect 

coefficient was approximated by shifting the indirect portion of the 

a vs. E profile of material #1 by The difference between the indirect 

bandgaps of material #1 and the alloy. Similarly, the direct 

coefficient was approximated by shifting the a vs, E profile of 

material #2 by the difference between the direct bandgaps of material #2 

and the alloy. This procedure only requires that the absorption 

coefficient vs, photon energy curves for both binary components be 

available, along with the bandgap vs, composition relationship (see 

Section 4,3,2) for the alloy. The results of this interpolation 

which is similar to that used by Hutchby [30], appear in Figures (5.5) 

and (5,9) in Chapter 5 for A1 Ga., As and GaAs., P , 
r x 1-x 1-x x 

Now, if the composition profile, C(X) is known, then, by the above 
procedure, a(A,x) is known, so that Equation (4,31) can be solved as 

N(A,x) = N(X,0) expE-/ X a(X,y)dy]„ (4.32) 

Thus, if each absorbed photon generates one electron-hole pair, the 
generation rate for a given wavelength can be written as 

G (A,x) = - — = cs(A,x) N(A,0) exp[-/ X a(X,y)dy], (4.33) 

Finally, the total generation rate for all values of X is 

G e (x) - £ T(A)N^(A)a(X ,x) exp£-/*a( A ,y )dy], (4.34) 

X 

Of course, if the semiconductor composition does not vary with 
depth, then Equation (4,34) reduces to the simpler form used in 
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calculating xhe generation rate profile for homogeneous solar cells 
[19]. 

Figure 4.3- shows the flow chart for the program --used . to calculate 
the 'generation rate profile for -analysis of variable composition cells „ 
Among the secondary outputs available as user options are: the 

dependence of antireflection layer efficiency on layer thickness; the. 
position dependence of the absorption coefficient; crystal surface 
and Si0 refractive indices; and 'the total optical current available as 
a function of ‘sample .thickness. The details of the numerical techniques 
used in the program can be found in Reference [19]. 







Figure 4.3 Flow chart for optical generation rate program. 
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5o PROPERTIES OF MATERIALS STUDIED 
5.1 Introduction 

At the present time, the computer program discussed in the preced- 
ing chapters has been used to analyse the performance of solar cells 

made of A1 Ga_ As- In Ga, As and GaAs, P - Since the results of 
x 1-x x 1— x 1-x x 

these analyses must be interpreted in light of the experimental data 
fed into the program 9 this chapter* presents the values used for the 
material parameters for these alloys- Every effort has been made to 
collect the latest, most reliable data in order to maximize the accuracy 
of the computer predictions. 

5.2 Material Properties of A1 Ga. As 

x 1-x 

One of the most promising alloys for solar cells is Al^Ga^^As. 
Since it has been used extensively in semiconductor laser research, its 
fabrication technology is well established and experimental data 'on its 
electrical and optical properties is relatively abundant. The lattice 
match between AlAs and GaAs is one of the best (less than .14% lattice _ 
mismatch), to be found among the material pairs suggested for solar 
cells so that recombination losses at any AlAs-GaAs heterojunction 
should be of minimum consequence E'29] - In addition, AlAs is an indirect 
semiconductor with a much larger bandgap (2.16eV) than GaAs (1-439 eV) 
which is a direct material- Such a combination permits the selection 
of a wide range of optical properties that should lead to interesting 
computer analysis results. 
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The following material parameters were used in the computer 
analysis of A1 Ga.. As solar cells . The subscript 1 refers to AlAs 
while the subscript 2 refers to GaAs„ 

The two sets ■ of relative dielectric constants used were [35]: 



All calculations were done for a temperature of 30O o K o The single 

trapping center -assumed by the Shoekley-Read-Hall recombination model 

in Equation (3,32) was set near the center of the bandgap regardless 

of alloy composition so that 

n.= p, = n. (5«1) 

1 1 i 

Of course the two parameters n^ and can easily be changed to fit 
the characteristics of any trapping level within the bandgap „ 

The procedure discussed in Sections (4.3,^.) and (4»3<,5i) for 
determining mobility, as- a function of composition and doping requires 



The following set of equations which represent ’curves fitted to 
available data (see Figures 5.1, 5,2 and 5; 3).' 

Holes,, data for ,GaAs used [373: 

2 

_ 380cm /V*sec 

^2.= V, = Cl+(3 . 17xl e- 17 cm 3 )N] 0 ” 266 

where N'= total impurity density 

Electrons, indirect valley, data for AlAs used [383: 


y . = f -i (N) = 
nl - nl 


165cm /V«s ec 


r , / 0 i - -.“17 3\» t *i0#13 

[It (8. 1x10 cm )N3 


Electrons, direct valley, data for GaAs used [373: 
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(5.2) 


(5.3) 


v n2 - f „2 <N> * 


7200cm /V-sec 


r 1 r r r -t 1n "17 3 \. T -l0 .233 

[If (5. 51x10 cm )N3 


(5.4) 


By curve fitting experimental data [31,33,34], the .following 
equations were used to, estimate the diffusion lengths of electrons and 
holes : 


L = 
n 


8yM 


-IQ Q 

ItN ( 8x10 /cm ) 


(5.5) 


3yM 


p ltN(1.2xl0~ 18 /cm 3 ) 
where N = total doping density. Then, by using Equations (5.2) and 
(5,4) to calculate the mobilities, the lifetimes required 'by the 
Shockley-Read-Hall model were estimated to be: 

.2 


(5.6) 


n 


no kT 


=5,3 nsec 


T yn 


po 


kT 

~T v p 


8.5 nsec 
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DOPING DENSITY (SITES/CM 2 ) 


Figure 5 «2 Empirical relationship between electron mobility and 
doping in AlAs as given by Equation (5,3) 0 







DOPING DENSITY (SITES/CM 2 ) 


Figure 5.3 Empirical relationship between electron mobility and 
doping in GaAs as given by Equation (5 . 4} „ 
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17 3 

when the total doping concentration was 4x10 /cm in the top n-type 
17 3 

layer and 2x10 /cm in the p-type base. 

The values used for the direct and indirect bandgaps of the two 
binary semiconductors were: 


AlAs: 

direct gap = 

E gdl 

= 2.95eV [39] 


indirect gap = 

E gil 

= 2.16eV [40] 


absolute gap = 

E gl 

= E . = 2.16 eV 

gil 

GaAs: 

direct gap = 

E gd2 

= 1.439eV [41] 


indirect gap = 

E gi2 

= 1.87eV [41] 


absolute gap = 

E g2 = 

: E ._ = l,439eV 
gd2 


The direct gap for the ternary alloy was calculated as a function of 
composition by [41]' 

E .(C) = 0.468C 2 + 1.042C + 1.439 (5.7) 

gd 

where C = mole fraction of AlAs. 

The indirect gap for the alloy was approximated linearly as 

E .(C) = 1.87 + 0.29C (5.8) 

gi- 
lt is apparent from the bandgap vs. composition curves of Figure 

5.4 that the direct-indirect transition occurs at about 44% AlAs, in 

approximate agreement with the transition point determined experimentally 

by others [42]. As discussed in Section 4.3.6 the electron affinity is 

assumed to depend on bandgap through a proportionality constant K . The 

cl 

values used for K and for the electron affinity in the two binary semi- 
conductors were [26,35]: 
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K = 0.85 
a 


A1AS! X -' 3,57eV o&^is' 


GaAs: x c2 = l+ * 07 eV 


OF PnnZ S 

P00ie QUALFIY! 


Specifically, these values -give electron affinity as a function of 
bandgap and thus composition by ‘Equations (5.7) and (5.8), as 


X = 4.07 - 0.85CE . 
c gd 


E 0 ) + E , - E 
g2 gd g 


(5.9) 


The values for the lattice constants required for determining interface 
state density as discussed in Section (4.4) are [29] 


o 

AlAs: a^ = 5.661A 

o 

GaAs: = 5.654A 

Surface recombination velocities at the two ends of the semiconductor 
crystal were chosen to be 

back surface: S = “cm/sec 

5 

front surface: S = 10 cm/sec 

(except for one computer run discussed in Chapter 6). 

The optical carrier generation rate G^ was calculated in the manner 
described in Section ”4.5 . Figure 5.5 shows the results of the inter- 
polation procedure described in that section for determining the absorp- 
tion coefficient as a function of both composition and photon energy. 

(or wavelength) . A similar procedure was used to estimate the refract 
tive index vs. - wavelength profile for an arbitrary alloy composition at 
the semiconductor - Si© interface. Briefly, the refractive index curve 
for GaAs was written 'as a function of photon energy, and then shifted 
by the difference between the ternary alloy bandgap and the bandgap of 
GaAs. The resulting curve for AlAs can be compared with the original 
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Figure 5 .6 
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curve for GaAs in Figure 5.6. This estimation procedure can be improved 
whed data concerning the refractive index of AlAs or A1 'Ga As becomes 

X 

available . 

All computer runs made as of the time of this writing have used the 
AMO solar irradiance data shown in Figure 2.2. By running the indepen- 
dent generation rate program with a range of SiO antireflection layer 
thicknesses, it was determined that maximum excess carrier generation 
occurs at a layer thickness of 0.07yM (±0.01pM) regardless of composition 
profile . 

The results of the computer analysis of A1 Ga As solar cells will 

X X"X 

be presented in Chapter 6. 

5.3 Material Properties of GaAs P 

±-'X X 

A second interesting semiconductor alloy for solar cell applications 
is GaAs P which, in many respects, resembles A1 Ga As. GaP, like 
AlAs is an indirect material having a bandgap considerably larger 
(2.262eV) than that of GaAs. Since the band structures of Ga,P and AlAs 
are similar, it is not surprising that their optical properties also 
resemble each other. Since GaP has a slightly larger absolute bandgap 
than AlAs , it might be expected than GaP would make a better window layer 
by absorbing a smaller fraction of the solar spectrum. However, it will 
be seen than the direct bandgap of GaP is considerably smaller than that 
of AlAs so that higher energy photons are absorbed at a greater rate 
than in AlAs. Given the spectral distribution of sunlight under AMO 
conditions the ultimate effect is greater attenuation in a GaP window 
than in an AlAs window of the same thickness. This fact will be shown 
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clearly by the, profile for the generation rate near the front of the 
cells discussed in Chapter 6. 

Probably the most serious problems inherent in GaAs P solar 
cells are associated with the fact that the lattice mismatch between GaP 
and GaAs (about 4%) is considerably greater than that of AlAs and GaAs 
C29j. This can be expected to lead to much more severe recombination 
losses at any GaP-GaAs hetero junction present in such a device. How- 
ever since one purpose of this research is ro clarify the importance of 
interface recombination to solar cell performance, this relatively poor 
lattice match makes GaAs.. P an interesting material for computer 
analysis. 

Although experimental data for the properties of GaAs P are not 

-L'-X X 

as readily available as data for A1 Ga. As , enough information has been 

X JL~X 

found to satisfy the requirements of the computer simulation program. 

In the following data, the subscript 1 refers to GaP, while the subscript 
2 refers to GaAs. Unless otherwise stated, the material properties per- 
taining to GaAs are the same as those used in the analysis of A1 Ga As 

X X""X 

listed in The previous section. 

The two relative dielectric constants for GaP were assigned values 
of C43H: 

low frequency! s = 11.1 
high frequency: e. 1 = 9.11 

The carrier effective masses for GaP in terms of electron rest mass 

m are: 
o 
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holes: m* = 0.772m [44] 

pl o 

electrons, direct valley: ra *c<fll = 0 *114111 £443 

indirect valley: m* . . = (0.34) (3) [45] 

J oil o 

Once -again, all computer runs were made for a temperature of 300°K and 
the Shockley-Read-Hall recombination parameters were given by Equation 
(5.1) for all composition values. The lifetimes used in the S-R-H model 
were once again: 

t = 5 . 3nsec , 
no 


t = 8.5nsec. 
po 

The doping dependence of mobility for electrons in the indirect 
valley used an empirical equation fitted to data for electrons in GaP 
[38]: 


y nl " f nl (N) 


200 cm, /V* sec 


r , _ ..-18 3. .,-,0.244 5 

11+ (1.5x10 cm )NJ 


(5.10) 


where, as before, N = total impurity density. The above equation is 
plotted for reference in Figure 5.7. 

The values for the direct and indirect bandgaps of GaP and GaAs 


were taken to be 


GaP: 

direct gap 

” 

E = 2.783eV 

gdl 

[46] 


indirect gap 

= 

E . = 2 . 262eV 

gll 

[47] 


absolute gap 

= 

E = E . = 2 . 262eV 

gl gil 


GaAs : 

direct gap 

= 

E = 1.425eV [46] 

gd2 



indirect gap 

= 

E ._ = 1.907eV [47] 

gi2 



absolute .gap 

= 

E . = E » 1. 425eV . 

g2 gd2 
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Although the electron , affinities for the two binary semiconductors are 
virtually the same, the composition dependence of bandgap for the 
ternary alloy means that -electron affinity cannot be taken .as 'Constant . 
From the discussion in Section 4.3.6 the dependence of electron 
affinity on . bandgap can be writxen explicitly as 

X = 4.07 - 0 . 384(E ,-E.-) + E ,-E . (5.13) 

c gd g2 gd g 

The laxrice constant for GaP was taken as [293 

= 5.451A , 

Once again the surface recombination velocities were* specified as: 

back surface: S = ra cm/sec, 

5 

front surface: S = 10 cm/sec (except for one run discussed 


in Chapter 6); 



MOLE FRACTION Gfi-P 


Figure 5, 


1 Bandgap of GaAs } as a function of composition [ u 6 9 47 ] 
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Finally, the optical generation rate for GaAs P was dexermined 

1 X X 

using the .same absorpxion coefficient and refractive index interpolation 
procedures used for A1 Ga As. The wavelengxh and conposition 
dependence of absorption coefficient and refractive index produced by 
these techniques can be seen in Figures (5. *9) and (5.10). Once again, 

AMO solar irradiance conditions were specified for all runs and an 
optimum SiO antireflection layer thickness of 0.07 yM was used. The 
computer predictions concerning the performance of GaAs P solar 
cells will be discussed in the next chapter. 

5.4 Material Properties of In Ga As 

X X“X 

Among the problems considered in this research is the determination 
of the relationship between substrate bandgap and solar cell efficiency. 
Although initial computer simulations for this purpose examined the 
performance of cells having substrates made of either GaAs, A1 Ga As 
or GaAs.. P , these materials have bandgaps greater than or equal to 

±~X X 

that of GaAs which is about 1.44 eV. In order to study the effects of 
smaller bandgaps , cells having substrates made of In Ga. _ As were 
subsequently analyzed. Unlike A1 Ga As and GaAs- P , In Ga As is 

X _L-X J.-X X X JL X 

a direct material for all values of x. The high absorption rate 

i 

resulting from these properties is well suited to a substrate, but not 
to a window layer. For this reason, the cells analyzed in this work 
having In Ga. As substrates have used AlAs windows. 

X X X 

One of the most . crucial parameters to be remembered during the 
interpretation of the computer simulation of AlAs-on-In Ga As solar 
cells is the lattice mismatch (about 7%) between InAs and both GaAs 
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and AlAs, which can be expected to produce significant recombination 
losses due to lattice mismatch in variable composition layers. The 
effects of bandgap changes and interface states introduced by alloying 
InAs with GaAs beneath an AlAs window will be discussed in the next 
chapter . 

Consistent with the notation established earlier, the subscript 1 
in the following data refers to InAs while the subscript 2 refers to 
GaAs. If not otherwise specified the parameters for GaAs are the same 
as those given in Section 5.2. 

The two relative dielectric constants for InAs were [35] : 
low frequency: e = 14.55, 

high frequency: = 11.8. 

The carrier effective masses for InAs were assigned values of: 

holes: m* = 0.41 m [35], 

pi o 

electrons, direct valley: m* = 0.024 m [35], 

J cdl v o 

indirect valley: m*.. = (0.39) (3)2^111 . 

* " cil o 

The effective mass for indirect valley electrons in InAs was set equal 

to that in GaAs for lack of further data and because In Ga. As is 

x 1-x 

entirely a direct material so that indirect valley electrons do not 
significnatly contribute to carrier activity. For the same reason, 
the doping dependence of mobility for indirect valley electrons was 
calculated from the empirical relationship used for such electrons in 
AlAs as given by Equation (5,3), As for other .materials analyzed in 
this work, the Shock ley-Read-Hall recombination parameters were given 
by Equation (5.1) and the lifetimes were specified as: 
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x =5.3 nsec, 
no * 

r = 8.5 nsec. 


po 


The values for the bandgaps of InAs are: 
direct gap = E - 0.359 eV 


indirect gap = E = 1-60 eV 


£ 353 , 

in* 


absolute gap = E,=E,, = 0.359 eV, 
^ gl gdl 


and the corresponding equations for the direct and indirect bandgaps of 
In Ga As are (see Figure 5.11): 

E gd (C) = '1.439-1. 41C + 0.33c 2 [7,35,413 (5.14) 

E .(C) = 1.87-0.27C (5.15) 

gl 

The parameters used to determine electron affinity as discussed in 
Section 4.3,6 are: 


K = 0,77 
a 

InAs:x c -^ = 4.90 eV [353, 

so that electron affinity was calculated from tl>e jbandgap as : 

X = 4.07-0.384(E ,-E _). (5.16) 

c gd g2 

The lattice constant for InAs is 


a = 6.058 A [293. 


Once again the optical generation rate for In Ga 1 As was calculated 
by the same interpolation procedure discussed at the end of Section 5.2. 



MOLE FRBCTION IN-flS 


Figure 5,11 Bandgap of In Ga As as a fancxion of comoov.ixion 

X X^X " 
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6. COMPUTER ANALYSIS RESULTS 
6.1 Introduction 

The computer analysis program has been designed te allow for the 
specification of a, wide variety of variable, composition solar cell 
structures. Device polarity can be either n-on-p or p-on-n with either 
an abrupt or a Gaussian doping profile. Doping levels can be 
arbitrarily specified, and, if desired, a heavily doped layer can be 
placed at the back contact to create a high-low junction for improved 
carrier collection and reduced dark current (see Reference [193). The 
composition profile can be arbitrarily defined by providing the position 
dependence of the relative concentrations of two binary semiconductors 
in a ternary alloy. Dark characteristics are obtained by selecting an, 
optical generation rate of zero, and characteristics in light can be 
found for either AMO or AM2 illumination conditions. Finally, an 
antireflection layer made of SiO can be assigned any-reasonable thickness. 

The simulation results discussed in this chapter include' Al^Ga-j _^As , 

In Ga, As and GaAs, P cells where the composition is either linearly 
x 1-x 1-x x 

graded or abruptly changed at a hetero junction and fs "constant in the bulk 
beneath the window layer. All cells analyzed up to the present time 
have used an 0.07 pM (near optimum) SiO antireflection layer and have 
been illuminated under AMO conditions. The structure and key performance 
parameters for twenty-nine of the more interesting devices can be found 
in Figure 6.1 and Tables 6.1 and 6.2. 

The first device shown in Table 6.1 is an n-on-p GaAs solar cell 
with no composition variation. The data obtained fer this "control" 




Figure 6.1 Structure common to all devices in Tables 6.1 and 







Table 6.1 Structure parameters and computer analysis results for severed. (Al,In) Ga As solar cells 
(see Fig. 6.1), under AMO irradiation. x -x 


pevice 

Alloy 

Comp. 

Het. 

Mole 

p-n 

Polarity 

S 

Inter- 

V ' 

j 

1 1 

Fill 

Peak 

# 


Change 

Depth 

Frac, 

Depth 


cm/sec 

face 

oc 

(V) 

SC 

2 

Factor 

Eff. 




(yM) 

GaAs , 
Bulk 

(yM) 



States 

(mA/cm ) 


(%) 



1 

GaAs 

none 

- 

1.00 

1.0 

2 

A1 J3a n As 
x 1-x 

linear 
(-to 35% 
AlAs) 

1.0 

1.00 

1.0 

3 

A1 Ga. As 
x 1-x 

linear 

1.0 

1.00 

1.0 

4 

A1 Ga. As 
x 1-x 

linear 

1.0 

1.00 

1.5 

5 

A1 Ga. As 
x 1-x 

abrupt 

1.0 

1.00 

1.0 

6 

A1 Ga, As 
x 1-x 

abrupt 

1.0 

1.00 

1.0 

7 

A1 Ga. As 
x 1-x 

abrupt 

1.0 

1.00 

1.3 

8 

A1 Ga. As 
x 1-x 

abrupt 

0.5 

1.00 

0.8 

' 9 

A1 Ga. As 
x 1-x 

abrupt 

0.1 

1.00 

0.4 

10 

Al Ga. As 
x 1-x 

abrupt 

0.1 

1.00 

0.4 

11 

Al Ga.. As 
x 1-x 

abrupt 

0.1 

0.95 

0.4 

12 

Al Ga. As 
x 1-x 

abrupt 

0.1 

0.90 • 

0.4 

13 

AlAs/In Ga.. As 
x 1-x 

abrupt 

0.1 

0.99 

0.4 

14 

AlAs/In Ga. As 
x 1-x 

abrupt 

0,1 

0.98 

0.4 

15 

AlAs/In Ga As 

X 1 -X 

abrupt 

0.1 

0.95 

0.4 

16 

Al Ga. As 
x 1-x 

abrupt 

0.1 

1.00 

0.1 

17 

Al Ga. As 
x 1-x 

abrupt 

0,1 

1.00 

1,4 


n-on-p 

10 

f™ 

.964 

17.8 ■ 

.822 

10.42 

n-on-p 

10 5 

No 

.999 

32.1 

.864 

20.46 


5 






n-on-p 

10 

No 

.999 

32.6 

.863 

20.75 

n-on-p 

10 5 

No 

,998 

32,4 

,863 

20.60 

n-on-p 

10 5 

No 

,998 

31,4 

.861 

19,95 

n-on-p 

io 5 

Yes 

,991 

31,4 

.862 

19.85 

n-on-p 

io 5 

Yes 

.998 

31.8 

.862 

20.21 

n-on-p 

io 5 

Yes 

.999 

32.3 

.863 

20,55 

n-©n-p 

io 5 

Yes 

.999 

32.8 

.867 

20,95 

n-on-p 

io 6 

Yes 

.999 

32.8 

.863 

20.86 

n-on-p 

io 5 

Yes 

1,050 

31.0 

.868 

20.83 

n-on-p 

10 5 

Yes 

1.103 

29,1 

.873 

20.62 

n-on-p 

10 5 

Yes 

.986 

33.3 

.864 

20,94 

n-on-p 

io 5 

Yes 

.973 

33.7 

.867 

20.92 

n-on-p 

io 5 

Yes 

.934 

34.8 

.865 

20.83 

p-on-n 

io 5 

Yes 

1.032 

29.8 

.850 

19.27 

p-on-n 

io 5 

Yes 

1,028 

32,7 

,857 

21,23 


Table 6.2 Structure .parameters and. computer analysis results for several GaAs P solar cells 
(see Fig. 6.1), under AMO irradiation. -x x 


Device 

# 

Alloy 

Comp. 
Change 1 

Graded . 
Width 
<«) 

Het. 

Depth 

<pM) 

Mole 
Frac. 
GaAs , 
Bulk 

p-n, 

Depth 

(pM) 

S 

cm/sec 

Inter- 

face 

States 

V 

oc 

(V) 

J 

sc 2 
(mA/ cm; 

Fill 

Factor 

Peak 
Eff . 
(-%) 

18 

GaAs. P 
ItX x 

abrupt 

- 

1.0 

1.00 

1.0 

io 5 

No 

.985 

30.64 

.862 

19.24 

19 

GaAs, P 

1-K X 

abrupt 

- 

1.0 

1.00 

1.0 

io 5 

Yes 

.831 

30,64 

J 799 

15.04 

20 

GaAs, P 
1-x X 

abrupt 

- 

1.0 

1.00 

1.04 

io 5 

c 

Yes 

.974 

30.59 

.859 

18.92 

21 

GaAs, P 
1-x X 

abrupt 

- 

.5 

1.00 

.54 

10 5 

Yes 

.975 

31.48 

.859 

19.49 

22 

GaAs, P 
1-x X 

abrupt 

- 

.1 

1.00 

.14 

10 5 

Yes 

.976 

32.43 

.859 

20.10 

23 

GaAs , P 
1-X X 

abrupt 

- 

.1 

1.00 

.14 

io 6 

Yes 

.961 

32.19 

'.857 

19.56 

24 ■ 

GaAs, P 
1-x X 

abrupt 

- 

.1 

0.95 

.14 

10 5 

Yes 

1.007 

30,28 

.870 

19.39 

25 

GaAs , P 
1-x X 

abrupt 

- 

.1 

0.90 

.14 

io 5 

Yes 

1.052 

28,28 

.870 

18.91 

26 

GaAs , P 
1-x X 

linear 

50.0 

.1 

1.00 

.14 

io 5 

Yes 

.971 

32.34 

.866 

20,21 

27 

GaAs , P 
1-x X 

linear 

200.0 

.1 

1.-00 

.14 

10 5 

Yes 

.981 

32,49 

.865 

20.40 

28 

GaAs, P 
1-x X 

linear 

600.0 

.1 

1.00 

.14 

io 5 

Yes 

.983 

32,51 

.864 

20.43 

29 

GaAs, P 
1-x X 

linear 1000.0 

.1 

1.00 

.14 

10 5 

Yes 

,983 

32,41 

.864 

20,37 
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cell can be used as a reference for comparison with the results of the 

simulation of variable composition .devices. Devices 2-10 are 

nAl Ga As-nGaAs-pGaAs with various p-n junction depths, window 

thicknesses and window composition profiles. The next five cells 

(11-15) have the same structure as device 9 except that the substrates 

are either ,A1 Ga, As (devices 11 and 12) or In Ga, As (devices 13-15) 

rather than ’GaAs alone. The last two cells in Table 6.1 represent a 

sample of the structures analyzed to determine the possible advantages 

of the p-on-n polarity. Although most of the present research .has dealt 

with .abrupt heterojunction devices, three graded A1 Ga 1 As cells 

(devices 2-4) and four graded GaAs P cells (devices 26-29) were 

simulated in order to estimate possible advantages of the built-in 

fields and broadened interface state distribution that accompany 

composition grading. As shown in the ( third column of Tables 6,1 and 

6.2, the grading was linear and graded layers ranged in thickness from 
o 

50A to 1 pM, 

Window thicknesses between O.J. pM and 1.0 pM were examined and 
numerous runs were made to. determine the optimum- placement of the p-n 
junction for .several values of window thickness. Although the p-n 
junction depth 'below the hetero junction was allowed to vary between 
0.0 pM and 2.0 pM, only the most .significant configurations appear in 
the tables. In an effort to determine the relationship between peak 
efficiency and substrate .composition, , several runs were made with a 
fixed, non-rzero mole fraction of either AlAs, InAs, or GaE in the bulk 
for comparison with similar structures -having substrates made entirely 

of GaAs. The results of the four most significant of these runs appear 

< 

in the tables for .devices 11-15, 24, and 25. 
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Since the major reason for using a wide bandgap window layer is to 

eliminate large surface recombination losses . it was not expected that 

the value of surface recombination velocity S would significantly affect 

performance in devices having such a window and so S was generally 

5 

assigned a value of 10 cm/sec. However, in order to verify these 

6 

expectations, two runs were made with a value of S ~ 10 cm/sec. The 
resulting performance data, shown in the tables for .devices 10 and 23, 
should be compared with devices 9 and 22, respectively, which differ 
from 10 and 23 only in surface recombination velocity . Finally, in 
order to measure the influence of interface states on performance, 
several cells were simulated (Devices 2-5 and 18) without allowance 
for interface recombination. The relationships indicated by xhe data 
in Tables 6.1 and 6.2 are discussed and illustrated in more detail in 
the following sections. 

6 .2 Interface States 

One of the foremost criteria used in selecting a ternary alloy 
system for variable composition solar cell use is the difference in the 
lattice constants of the two binary constituents. As discussed in 
previous chapters, the closer the lattice match, the fewer xhe interface 
states ax any heterojuncxion and the lower the density of dangling bonds 
in any graded layer created by spatially varying the alloy composition <> 
The basic reason for using a wide bandgap window layer is to reduce 
carrier recombination losses at surface imperfections. If the interface 
between the wide gap window and the smaller gap substrate is 
characterized by a high interface recombination velocity due to severe 
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lattice mismatch, surface recomb ination losses may simply be “replaced by, 

comparable interface recombination losses. 

The small lattice mismatch (about 0.14%) between AlAs and GaAs makes 

The A1 Ga As system a -prime candidate for heteroj unction solar cell use. 
X * 

4 

Interface recombination velocities on the order of 10 cm/s -have been 

reported [28] for abrupt AlAs-GaAs hetero junctions while reasonable 

values for the surface recombination velocity of GaAs solar cells 

6 7 

generally range from 10 cm/sec to 10 cm/sec [31,32]. By using Equation 

“8 

(4.25) and taking the lattice constants of AlAs and, GaAs as 5.661x10 cm 

and 5.654x10 cm, respectively [29], the density of dangling bonds at an 

AlAs-GaAs heterojunction with a [100] orientation can be estimated to be 
12 2 

about 3.1x10 /cm . 

The more severe lattice mismatch (about 4%) between GaP and GaAs, 
on the other hand, means that the GaAs P system can be 'expected to 

< 

suffer more 'serious interface recombination losses than the ,A1 Ga, As 

x 1-x 

system, especially in abrupt hetero junction structures. By calculations 
similar to those used for A1 Ga, As, and taking the lattice constant of 

X _L X 1 

— 8 

GaP as 5.451x10 cm, the surface density of interface states at a GaP- 

13 2 

GaAs heterojunction should be about 9.5x10 /cm , or about 30 times , the 

density at, an AlAs-GaAs hetero junction. If the capture coefficients 

(see Section 4,4) in A1 Ga, As and GaAs, P are similar, then an inter- 

x 1-x 1-x x 5 

5 

face recombination velocity of about 5x10 cm/sec can be expected at 'a 
GaP-GaAs he,terojunction. Since such a value is comparable to surface 
recombination velocities found in GaAs cells, the prospects for building 
efficient GaAs P solar cells would appear to be relatively poor. 

J-~X X 
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However, computer analysis of such cells can be useful in several 
ways. First, it can aid in assessing the relative importance of inter- 
face states in controlling peak efficiency by permitting a comparison 
with similar analyses of cells made of alloy systems with better inter- 
face properties. Second, it can add to the body of evidence needed to 
determine whether the GaAs P system is as poorly adapted to solar 
cell use as first order considerations seem to indicate. Finally, com- 
puter simulation can test the ability of such techniques as composi- 
tion grading and separation of the p-n .junction from the heferojunction 
in reducing interface recombination losses. 

The effect of interface states on AlAs-GaAs hetero junction cells 
can be seen by comparing devices 5 and 6 .which are identical except 
that the analysis of device 5 omitted the interface states. Clearly, 
the . degradation of peak efficiency (19,95,% vs, 19,85%) caused by inter 
face recombination cannot be considered a major factor in controlling, 
the performance of A1 Ga As cells. The effect should be even less 
significant in graded composition cells, and in cells having the p-n 
junction placed slightly beneath the hetero junction. These considera- 
tions will be discussed in later sections of this chapter. 

The effect of interface states on abrupt GaP-GaAs cells was, as 
expected, found to be much more pronounced. Device 18 represents a 
structure with both the heterojunction and the p-n junction placed 1 pM 
beneath the front surface of the crystal, but with no interface states. 
Device 19, on the other hand, shows the results of analyzing the .same 
structure by including the modeling of interface recombination as 
discussed in Section 4,4, The reduction in peak efficiency from 



105 


19,25% to 15 o04% results from the large increase in dark current 'caused 
by recombination of carriers within the depletion region. ' As Figure 6,2 
shows, the introduction of interface states (device 19) increases dark 
current by more than two orders of magnitude when the heterojunction 
coincides, with the ,p~n junction. As will be discussed in the next 
section, this increase in dark current can be significantly reduced 
by placement of the p-n junction slightly below the he ter ©junction. 

The effect of an increase in dark current is a greater loss ‘ in 
maximum power delivered to any external load as more internal forward 
bias current flows for any given terminal voltage. The softening of 
the knee of the J-V characteristics reduces the curve factor (from 
0,862 to 0,799 in devices 18 and 19, respectively). In addition, the 
increased dark current reduces the open circuit voltage , which can be . 
defined as the voltage at which internal bias current equals the current 
collected from optically generated .carriers . 

It is clear that interface recombination can seriously degrade the 
performance of heterojunction cells made of semiconductor materials that 
would otherwise be well suited to solar cell use. The next two sections 
discuss structure modifications that can help minimize the .losses 
associated with interface states. 

6.3 Optimization of p-n Junction Placement 

The determination ©f the best location for the p-n junction in a 
solar cell is not 'a simple problem, especially for variable composition 
structures. Of course the p-n junction should be placed where the .built- 
in field around it can most .effectively collect optically generated 


CURRENT DENSITY CAMP/CM 2 ) 
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DEVICE 20 
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Figure 6„2 Dark currant density vs„ vcdtage for devices 18, 19, and 

20 (Table 5 „ 2) showing fh.- of fact of interface recombination 
and p-n junction displ^comant. 
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carriers before they recombine ax the surface, at the hetero junction, 
or in the bulk,, As mentioned earlier, a , high concentration of interface 
states within .the junction depletion region can increase dark current 
and significantly reduce peak efficiency. 

As a first order approximation it is often assumed that all carriers 

generated within one diffusion length of the p-n junction are collected 

by the built-in field 0 By applying Equations (5.5) -and (5.S) ~to the 

17 3 

devices 'having an n- type -layer of doping density 4x10 /cm on top of a 

17 3 

p-type substrate of doping density 2x10, /cm (see Figure 6.1), the 
minority carrier diffusion lengths can be estimated at about 7 yM for 
elecxrons and about 2 yM for holes. In other words, the effective 
collection width around the ,p-n junction in n-on-p cells is approximately 
2 yM toward the front surface (in n-type material) and about 7 yM into 
xhe substrate (p-type material) . For cells wixh the opposite polarity, 
but the same doping concentrations (see Figure 6.1 and devices 16 and 17 
in Table 6.1), xhe minorixy carrier diffusion lengths are abeux 2.5 yM 
for holes and 6 yM for electrons, so that the ( effective collection region 
expends approximately 6 yM above the p-n junction (in p-type material) 
and about 2.5 yM into the substrate (n-type material) . Since the 
effective .collection region above xhe p-n junction is about three times 
greater in p-on-n cells than in n-on-p cells (for the same doping levels 
in each layer), it can be expected that the optimum junction placement 
will be deeper for the former than for the latter polarity. 

However, since xhe doping change at xhe junction for either polarity 
is abrupt, xhe depletion region is -highly localized and is typically, 
about 0.1 yM wide. The region of peak generation for an abrupt 
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hererojuncrion cell is also localized to the region just beneath the 
hetero junction, In fact* computer calculations for abrupt AlAs-on-GaAs 
and GaP-on-GaAs cells indicare that over 80% of -the total .optical 
carrier generation occurs in the rop 1 yM of rhe GaAs substrate (see 
Figures 2.9 and 6*3), Thus it can be expected that the optimum location 
of the p-n junction may be much closer to the heterojunction than 
indicated by the diffusion lengths of minority carriers alone. 

The results of computer runs made expressly to .determine the best 
location of the p-n junction are shown in Figures 6.4, 6.5 and 6.6. 

The performance parameters for the seven structures 'that correspond 
to the peaks of rhe curves in These figures are lisTed in. Tables 6.1 
and 6,2 as devices 7,8.9 ,17, 20 ,21 and 22. Figure 6.4 shows The 
dependence of efficiency on p-n junction depth into the GaAs substrates 
below an n-type AlAs window of various thicknesses. For all Three 
values of window thickness shown. The opTimum placement; of The p-n 
juncrion occurred at approximately 0.3 yM beneath the abrupt hetero- 
junction. For more shallow placements the increase in dark current due 
to interface states with the depletion region caused a slight drop of 
about 0.3 percentage points in peak efficiency. At greater depths the 
collection efficiency of the junction suffered from the separation of 
the highly localized regions of carrier generation and built-in field. 

Figure 6.5 indicates, as expected from diffusion length .estimates , 
that the best location for the p-n junction in pAlAs-p'GaAs-nGaAs cells 
is deeper than in similar n-on-p structures. For a 0.1 yM AlAs window 
and the doping levels of Figure 6.1, a peak efficiency of 21.23% was 
obtained by placing the junction 1.3 yM below The window (device 17) 
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Dcpcnder^'i 3f pedl efficiency on p-n junction depth below 
the hererojunct izn for abrupt nALAs-nGaAs-pGaAs cells <, 
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Figure; 8„6 Dependence. of peak efficiency on p-n junction depth below 
the her*r o junction fob abrupt nGaP-nGaAs-pGaAs cells. 
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compared with an efficiency ©f 20.95% in the corresponding n-on-p cell 
(device 9) with the junction optimized at a depth of 0.3 yM beneath the 
window. Since the open circuit voltage of device 17 is slightly greater 
than that of device 9 while the short circuit currents are almost equal, 
the higher efficiency of the p-on-n cell is due to a lower dark current 
rather than a higher collection efficiency. Although it, appears that 
p-on-n solar cells may have slightly higher .efficiencies than cells of 
the opposite polarity, the difference is so small that the choice will 
more likely depend on other, factors such as process technology, and 
radiation tolerance. 

The relationship found for abrupt GaP-on-GaAs n-on-p cells, shown 
in Figure 6.6, shows more variation than that for AlAs-on-GaAs cells. 

The best location for the p-n junction was found to be 0.04 pM beneath 
the wide bandgap window. From computer simulation data, the distance 
between the p-n junction and the edge of the depletion region on the 
n-type side was found to be close to 0.04 pM, so that maximum efficiency 
can be obtained by placing the heterojunction (with its high concentra- 
tion of recombination centers) near the edge of the depletion region. 

If the p-n junction is placed closer to the window layer, peak effi- 
ciency drops rapidly until a value of only about 15% was calculated 
when the het ere junction and p-n junction coincide. Once again, this 
drop is due to an increase in dark current as shown by the dark J-V , 
characteristics for devices 19 and 20 in Figure 6.2. 

From this -analysis it appears that by careful positioning of the 
p-n junction, GaP-on-GaAs heter ©junction solar cells can be made to 
. perform almost as efficiently as similar AlAs-on-GaAs cells in spite 
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of the inferior lattice match of the GaP-GaAs system,, However, such 
accuracy naturally calls for expensive fabrication techniques. In fact, 
present technology would be severely strained t© reproducibly make 
devices corresponding to the peak region of Figure 6.6. It seems more 
reasonable instead to investigate the possibility of controlling inter- 
face state losses by spreading the composition change over a finite 
distance as suggested in the next section. 

6.4 Composition Grading Effects 
Theoretical calculations by other workers t 30 ,49 3 have indicated 
th^t built-in fields due to bandgap grading in variable composition 
solar cells can improve minority carrier- collection efficiency by 
reducing the surface recombination losses. For example. Figure 6.7 
shows the equilibrium band diagram for device 3 which is an n-on-p 
cell with a wide bandgap window, linearly graded from AlAs at the front 
to GaAs at the p-rn -junction 1 pM below the surface. According to 
Equation (2.8) such a composition profile should create an effective 
fiqld for holes in the n-type graded layer that will sweep them away . 
fpqm the front surface and toward the p-n junction where they can be 
collected. Of course, the importance of such a process depends on the 
density of excess carrier generation in the window,' which in turn is 
essentially controlled by the bandgap. For this reason, it is not 
immediately clear whether a constant wide bandgap window of AlAs is 
preferable to a graded window of A1 Ga As. In other words, it is 
pot obvious that the built-in field is effective enough to compensate 
for the fact that more minority carriers are generated farther from the 
junction in a graded layer than in a constant wide bandgap window of AlAs. 
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By comparing devices 3 and 5, which are identical, except for the 
composition profile within the window, it is apparent that the gain in 
peak efficiency obtained by using a graded window is slightly less than 
1 percentage point (20,75% vs, 19,95%). Although such an increase is 
not extraordinary, it is possible that significantly greater improve- 
ments can be made 'by examination of other, possibly nonlinear, grading 
profiles . 

In GaAs P devices, the relationship between composition grading 
and peak efficiency is complicated by the poor lattice match of the 
GaAs-GaP system. In order to 'determine this relaxionship , numerous 
simulations were made for n-on-p cells with the p-n junction placed 
at 0.14 uM beneath the surface and the window layer linearly graded 
over a variable width with grading terminated in GaAs at a depth of 
0,1 pM (see Figure 6.1). The results of four such simulations appear 
in Table 6.2 for devices 26-29 and the dependence of peak efficiency 
on graded layer width, is plotted in Figure 6.8. The initial gain in 
efficiency, as the graded width is increased from zero, is probably due 
to the widened distribution of interface states which had previously 
been concentrated at an , abrupt hetero junction. Such a distribution 
increases carrier lifetime in the region of changing composition result- 
ing in greater collection efficiency. Of course, the built-in fiefd due 
to bandgap grading also improves the collection efficiency, just as it 

did for A1 Ga. As cells discussed earlier, 
x 1-x 

However, as the graded width is increased beyond 0.04 jiM,.peak 
efficiency falls as more carriers are generated further from the p-n 
junction and progressively fewer holes survive the short lifetimes 
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that characterize the graded, GaAs^ region. Thus, it appears that 1 
grading linearly from GaP to GaAs over about 0,04 pM can, yield an 
improvement of about 0.4 percentage points (20.5% vs. 20.1%) in 
efficiency when compared with an abrupt GaP-GaAs composition profile 
when the ,p-n junction is placed at the -optimum depth for an abrupt 
composition profile. 

As discussed in the previous section, the, placement of the p-n 
junction is very critical to the performance of GaAs.. P solar cells. 

1 X-X X 

Therefore it would not be surprising to find that further gains could 
£e obtained by optimizing the location of the p-n junction, in a . 

GaAs. P cell in which composition is linearly graded over about 

X— X 4 

0.Q4 pM, Since the graded profiles 1 of devices ,26-29 shift carrier 
generation slightly toward the .surface , it would seem reasonable to 
expect the optimum p-n junction depth, to be slightly closer to the 
surface as well. Some, idea. of the possible increase in efficiency 
associated with optimization of the p-?n junction placement can be 
obtained from examination of Figure '6.6, which shows that a change of 
only 0.1 pM can result in a change of 1 percentage point in peak 
efficiency, even when the depletion region does not overlap the variable 
composition region. Once -again, however, such accuracy would be diffi- 
cult to achieve under the limitations of present technology. 

6.5 . Effect of Window Thickness on Efficiency 
The use of a wide bandgap window layer affects solar cell 
performance in two important ways. > First, it separates surface 
imperfections from the regions of high carrier generation rate just 
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beneath the window. The improved carrier collection efficiency that 
results ,frotn This separation is -the main advantage of the window layer. 
Second, since no practical .wide bandgap material is completely trans- 
parent to the enxire solar spectrum, the absorption of some photons, 
essentially those toward the .shorter wavelength end of the specxrum, 
is unavoidable. The thicker the window, the greater the absorption 
and the smaller the efficiency of collection of electron -hole pairs 
generated within the window. In other words, a trade-off exists 
between reduction of surface recombination losses and an increase in 
ordinary carrier recombination within the window layer itself. 

In an effort to determine the optimum window' layer thickness, 
several computer analysis runs were made for abrupt heterojunction 
AlAs-on-GaAs and GaP-on-GaAs solar cell structures. As discussed 
in Section 6.3, the depth of the p-n junction was varied in order to 
establish the best location at each of three wipdow thicknesses: 

O.lyM, 0.5uM, and l.OpM. Figure 6.9 shows the relationship between 
peak efficiency and window thickness for abrupt hetero junction AlAs- 
on-GaAs solar cells in which the placement of the p-n junction has 
been optimized. This relationship is similar to that reported by Hovel 
and Woodall in 1973 [583 from their comparer analysis of similar 
A1 Ga. As structures. Calculations by Hovel [54j also led to the same 

X J.™X 

general .conclusions. Although the computer analysis indicates that 
maximum efficiency occurs for window thicknesses less than Q.lpM, 
practical fabrication problems associated with thinner windows may 
outweigh any gains in efficiency that could be achieved. 
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Dependence of peak efficiency on heterojunction depth 
for abrupt AlAs-GaAs and GaP-GaAs solar cells, with 
optimum placemanx of The p~n junction. 
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Because of the presence of surface states,' as modeled by the surface 

recombination velocity S, peak efficiency can be expected to fall as the 

window thickness approaches zero. Since the surface recombination 
5 

velocity of 10 cm/sec specified for these computer simulations corre- 
sponds to relatively high quality surface characteristics, the-results 
might appear to be somewhat optimistic. In order to test the effect of 

a higher surface recombination velocity, a computer analysis was made - 

0 

for both alloy systems specifying a value of 10 cm/ sec -for S while using 
the same thin window structure used in devices 9 and 22. The results 
shown in Tables 6.1 and 6. '2 -as devices 10 and 23, indicate that such a 
ten-fold increase in S should reduce the .peak efficiency of, an AlAs-GaAs 
cell from 20.91% to 20.86%, and reduce the peak efficiency of a GaP-GaAs 
cell from 20,10% to 19.56%, These small reductions 'are due to a reduced 
carrier collection efficiency leading to slightly lower values of short 
circuit current. However, it is apparent that even a. very -thin window 

k 

' t ' "* ' 

layer (O.lpM) effectively eliminates surface recombination as a dominant 
loss mechanism for reasonable .values 'of S. 

6.6 Bu-lk Bandgap Opximization 

Although it is well, known that the optical .properties of GaAs are 
better suited to solar energy -conversion than tho.se of silicon [4,21, 
23,54] (see Figure 6.10)- it is not clear that GaAs is the ;best ma.terial 
available, especially when the full- range-, of electrical and optical 
characteristics are examined in a complete solar cell ’Structure under 
AMG conditions . Specifically, the bandgap (l.44eV) of a GaAs ’ substrate 
may ‘be loo small or too large to produce the highest peak efficiency. 
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Figure 6.10 Dependence of conversion efficiency limit on bandgap 
for AMO and AM2 irradiance conditions [54], 
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In general, the wider the substrate bandgap, the higher the open circuit 

voltage and the smaller The dark current. But as the bandgap is widened, 

the substrate is made more transparent to a larger share of the solar 

spectrum so that fewer carriers are generated toward the front of the 

cell near the p-n junction, resulting in a lower short circuit current. 

In other words, adjustment of the substrate bandgap involves a trade-off 

between open circuit voltage and short circuit current. 

The computer simulation of devices 9 and 11-15 in Table 6.1 and 

devices 22, 24 and 25 in Table 6.2 were made to determine the effect on 

efficiency of adjusting the composition of the substrate in abrupt 

he ter o junction (Al,In) Ga. As and GaAs.. P solar cells. In cells 11 

X -L“X X 

and 12, a 0.1 pH AlAs window covers an A1 Ga, As substrate of uniform 
5 x 1-x 

composition having a larger bandgap than GaAs, A 0.1 pM AlAs window 
is also used in devices 13-15, but the substrate i consists of In Ga. As 
which has a smaller bandgap than GaAs. The peak efficiencies for these 
cells are plotted as the solid curve in Figure 6.11. It is seen that a 
GaAs substrate is superior to one made of either A1 Ga As or 
In Ga. As. However, the inferiority of In Ga. As is due to the severe 
lattice mismatch (about 7%) between InAs and GaAs rather than the 
reduction in bandgap. When -the lattice constant of InAs was assigned 
a value equal to that of AlAs, the computer analysis predicted the 
performance indicated by the dashed curve segment in Figure 6.11, which 
peaks at a bandgap of about 1.41 eV. In other words, a substrate of 
In AO Ga OQ As would perform slightly better than a GaAs substrate if it 

o 1/4 o 3 0 

were not for the poor lattice match. However, it should be noted that 
the p-n junction placement was optimized only for the structure of 
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Figure 6.11 Peak efficiency vs. substrate composition in abrupt 
hetero junction AlAs-on-Al Ga As and AlAs-on-In Ga 
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cells with an O.lyM window and a p-n junction depth of 
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device 9 by varying the junction depth while holding other structure 
specifications constant- At the present time, this procedure has not 
been done for- any structure having a percentage of AlAs or InAs in the 
substrate- More specifically, because the introduction of AlAs into 
the substrate widens the bandgap and therefore reduces the absorption 
coefficient for photons near the band edge, it should be expected that 
the best p-n junction placement should be somewhat deeper for (Al,Ga)As 
substrates than for a GaAs substrate. Conversely, an (In,Ga)As sub- 
strate probably calls for a shallower p-n junction. Since the 
generation rate distribution is very sensitive to material bandgap 
changes and the optimum p-n junction location in turn depends on the 
generation rate distribution, the trend shown in Figure 6.11 does not 
prove that further small gains in efficiency cannot be obtained by 
alloying a GaAs substrate with a small percentage of AlAs, or more 
likely, InAs. In any case, the potential improvement is so small that 
these computer simulations essentially reinforce earlier estimates that 
the bandgap of GaAs is near optimum for AMO conditions C543« 

Similar results were obtained from the simulation of cells having 
GaP alloyed ..with GaAs in the substrate beneath a 0.1 pM GaP window. 

The data for devices 22, 24 and 25 indicate that when the placement 
of the p-n. junction is optimized for an abrupt GaP-on-GaAs cell 
(device 22) the subsequent introduction of 5% GaP (increasing the band- 
gap from 1.44 eV to 1.48 eV) in the substrate reduces peak efficiency 
from 20.10% to 19.39%. When the GaP level is increased to 10% 

(E = 1.54 eV), the peak efficiency falls to 18.91%. Although this 
§ 

trend, shown in Figure 6.12, does not encourage the use of a substrate 
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Figure 6,12 Dependence of peak efficiency <?n substrate composition 
in abrupt heterojunction GaP-GaAs.. P cells with 
hetero junction depth of O.lyM and p-n junction depth 
of 0.14yM„ 





ORIGINAL PAGE If 
0F p OOR QUALITY 


127 


with a bandgap wider than GaAs, it must again be remembered that the 
p-n junction in devices and 25 may nox be placed deep enough for 
maximum efficiency. The stronger dependence found for Ga(As,P) 
substrates compared with (Al,Ga)As substrates may be due to the stronger 
relationship between p-n junction depth and efficiency found for the 
GaAsP system in Section 6.3 (compare Figures 6.4 and 6.6). 

If the performance of (Al,Ga)As, (In,Ga)As or Ga(As,P) substrates 
without optimum p-n junction placement had been found clearly superior 
to optimized devices with GaAs substrates, the analysis predictions for 
alloy substrates would have been favorable. However, in view of the 
results that were obtained, there seems little room for improvement of 
the GaAs substrate by alloying it with other materials to adjust the 
bandgap . 


6.7 Interpretation of Analysis Results 
It is important to view the preceding analysis results with 
consideration given to the techniques that produced them and the fabri- 
cation technology that is currently available to create the structures 
they call for. First of all, some of the material parameters, in 
particular excess carrier lifetimes, represent the characteristics of 
relatively high quality material. Analyses using such values may lead 
to slightly lower dark currents, and thus slightly higher open circuit 
voltages, fill factors and efficiencies than have been reported in 
experimental studies. Secondly, in comparing computer predictions with 
experimental data, it should be noted that most practical variable 
composition cells tested up to the present time have used less than 100% 
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AlAs or GaP at the surface. However, because early computer simulation 
data indicated that 100% AlAs or GaP at the front produced the highest 
efficiencies, most of the structures reported in this work specify such 
surface compositions. It is believed that the technical problems that 
discourage the use of AlAs and GaP as surface materials will soon be 
overcome . 

The difficulty of obtaining very thin window layers (less than 
O.lyM) may discourage the implementation of the optimum structures 
indicated in Section 6.5. However, the basic purpose of computer simu- 
lation is to permit reasonable estimation of the potential benefits of 
various solar cell structures so that costs can be weighed against 
these benefits. 

Finally, sheet resistance effects, which depend on contact 
geometry and are not incorporated into the analysis program, may become 
significant when the p-n junction is placed near the surface. The 
performance parameters reported in this work are for a one dimensional 
cell uncorrected for contact coverage and surface layer sheet resistance 
losses because all current has been assumed to be perpendicular to the 
front surface. In order to avoid the prohibitive optical losses 
associated with continuous contact coverage, practical cells generally 
use finger contact geometry which requires some lateral current flow 
near the surface. 

Although computer methods to correct simulation results for contact 
geometry are available [19], they were not considered necessary to the 
present research because the basic aim was to assess the relative 
influence of structure and material parameters on efficiency. Practical 
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contact areas Typically account for a 5-10% reduction in short circuit 
current and efficiency while sheet resistance may -reduce efficiency by 
about 0.5 percentage points. In other words, for a complete solar cell, 
the efficiency values reporred in earlier sections must be reduced by- 
about 1.5 percentage points so that a 20% cell should be considered an 
18.5% cell -after contacrs have been included in the estimate. 
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7. SUMMARY AND RECOMMENDATIONS 
FOR FURTHER RESEARCH 

7.1 Summary 

The aim of the present research was to develop a computer program 

capable of analyzing hetero junction and graded composition solar cells 

and to apply that program to a number of proposed cell structures using 

feasible semiconductor materials. The. program , now complete, accepts 

a wide variety of device structures and, by numerically solving a set of 

differential equations, is capable of simulating the behavior of .solar 

cells made of a ternary alloy between two binary semiconductors. Up to 

the present time, the program has been used to analyze variable 

composition solar cells made of A1 Ga, As, In Ga, As and GaAs, P . 

x 1-x x 1-x 1-x x 

The highest efficiencies calculared for A1 Ga As devices were 21.23% 

X X X 

for a pAlAs-pGaAs-nGaAs cell(device 17 in Table 6.1) and 20.95% for an 
nAlAs-nGaAs-pGaAs cell (device 9 in Table 6.1). When the values for 
peak efficiency obtained from these computer simulations are adjusted 
for contact geometry losses (typically a reduction of about 1.5 
percentage points) they are in good agreement with earlier less exact 
calculations [30,49 ,54, 58j. These calculated efficiencies are somewhat 
higher than some values reported for experimental variable composition 
cells such as 'the 14.7% devices built by Hovel and Woodall in 1973 [12 
This is not unexpected because the material parameters used in the 
present computer studies characterize relatively high quality semi- 
conductor material rather than the typical properties of most 
contemporary experimental cells. Also, the computer simulations 
reported in this work generally specified 100% AlAs or GaP at the surface, 
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while most experimental cells presently retain seme percentage of GaAs 
throughout. These calculations have been intended to determine the 
potential performance .levels of A1 Ga. As,' In Ga As and GaAs. P 
solar cells rather than explain individual experimental test ‘results. 
However, the latest results available at the time of this writing 
Include a report of an A1 Ga As cell with a measured AMO efficiency 
of 18.5% including contact ■ losses [61]], which is less than a percentage 
point below the best efficiencies calculated for such devices in the 
present work. 

Several relationships between device sxructure and performance 
have been, indicated by the analyses made up to the present time. The 
effect of heterojunction interface states was examined by comparison 
of computer , simulations with and without interface recombination. The 
resulting losses at abrupt AlAs-GaAs interfaces were minor. However, 
as expected from the severe 'latxice mismatch., the peak efficiency of 
GaP-GaAs heterojunction cells suffered far more from interface recom- 
binations which reduced efficiency as much as 5 percentage points when 
the p-n junction coincided, with the hetero junction. Subsequent computer 
simulations indicated than these losses could be largely offset by 
moving the p-n junction deeper so that the GaP-GaAs interface no longer 
occurred at the center of the depletion region. A further increase in 
efficiency of about 0.5 percentage points for GaAs^_ x P x cells .was 
obtained by linearly grading the change from GaP to GaAs over a distance 
of about 0.04yM. Grading over greater widths- led to significant 
recombination losses within the window layer. Similar, but, far less 
pronounced effects were observed in A1 Ga. As .simulations. It -was 

X 1™X 



132 


found "that p-on-n A1 Ga, As cells had smaller dark currents than 

similar n-on-p devices, allowing them re achieve slightly higher 

efficiencies and open circuit voltages. 

The success ©f wide bandgap window layers in reducing surfade 

recombination losses was corroborated by computer simulation of devices 

that differed only in surface recombination velocity. An increase in S 
5 6 

from 10 cm/sec to 10 cm/sec reduced the peak efficiency of AlAs-GaAs hells 
by less than 0.1 percentage points and thax of GaP-GaAs cells by less than 
0.6 percentage points. Since these results were obtained for window 
layers O.lyM thick, it is apparent that even very thin windows can 
virtually eliminate surface recombination as a dominant performance 
factor. 

Finally, computer simulation of devices having A1 Ga. As, 

X i.*”X 

In Ga As and GaAs.. P (0<x<0.1) substrates indicated that significant 

X i"X 1“X X 

gains in performance could .not be expected from use of substrates with 
wider or smaller bandgaps than GaAs . However , since the placement of 
the ,p-n junction was not optimized for such structures, it would not -be 
inconsistent with these calculations to find that small increases in 
peak efficiency could be obtained by careful positioning of the p-n 
junction, especially in the ,case of A1 Ga.. As and In Ga.. As substrates. 
In ’addition, the complex dependence of efficiency on substrate properties 
means that the use of materials other than GaAs cannot be ruled out 
simply on the basis of computer simulation of these alloy systems. 

The computer analysis program has giyen reasonable results for both 
abrupt and graded’ composition structures, and has been useful .in 
estimating the potential advantages of various configurations , 
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The detailed data which the program produces has been helpful in 
assessing the relative importance of the numerous inx enacting phenomena 
that control solar cell performance 0 

7.2 Recommendations for Further Research 

Although the computer analysis program is now essentially complete, 
its full potential has -not yet been utilized. Its capabilities permit 
the examination of several design options which have not been considered 
in the present work. Among the options which would require virtually no 
modification of the program are: 

1) Use of a Gaussian Cor other) doping profile 

2) Comparison of different doping levels 

3) Further .examination of p-on-n solar cells 

4) Optimization of p-n junction placement for all devices analyzed 

5) Application of irradiance conditions other than AMO, in 
particular AM2 conditions for terrestrial applications 

6) Use of antireflection layers other than SiO (for example TiO ) 

X 

7) Consideration of nonlinear composition grading to optimize the 
generation rate profile for improved collection efficiency 

8) Analysis of devices made of other alloys such as GaAs^ Sb 

-L X X 

9) Determine the effects of lifetimes shorter than those , examined 
in this work, i.e. consider use of lower quality, cheaper 
material. 

The study of several 'Other possibilities would require some minor 
program modifications, but could yield useful results. Among these are: 

10) Improved lifetime modeling, to relate lifetime di-rectly to 
doping levels and mobilities 

11) Design of a numerical step distribution within the analysis 
program that would relate positional step size solely to the 
slope of the three independent variables :Ji, <j> , and <p . 

(see Section 3.4, and References 13 and 19) n ^ 
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12) Performance under multi-sun and high temperature conditions 
such as those found in concentrator solar cell systems. 

Of course, the .input data concerning material parameters should 
continually be updated to maintain consistency with the latest available 
experimental figures. Finally, it is hoped that these computer studies 
will be followed by attempts to achieve the predicted performance levels 
in actual - working devices. After all, computer simulation is intended 
to provide guidance for experimental research. 
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